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PART I. 

Introduction. 

1. If in a gravitating body there occurs a displacement which involves alteration of 
density, there must be a tendency for the material to move towards the places where 
the density is increased, and away from the places where the density is diminished. 
The effect of this tendency, if it were not held in check, would be to accentuate local 
vol, ccvn. — a 417, z 2 31.5.07 
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alterations of density. In any body the tendency is partially held in cheek by the 
elasticity of the body, and, in particular, by the elastic resistance which the body 
offers to compression. If this resistance is sufficiently great, the body is stable, in 
spite of the tendency to instability which arises from gravitation. It is important to 
determine the conditions of stability for bodies of various forms and constitutions, 
with various distributions of density. The problem of the stability of spherically 
symmetrical configurations of a quantity of gravitating gas has been investigated by 
J, H. Jeans,* and he has drawn from his investigations some interesting conclusions 
in regard to the course of evolution of stellar and planetary systems. In a subse- 
quent memoir t he proceeded to investigate a similar problem in regard to gravitating 
bodies of a more coherent character. A gravitating solid body, such as a planet may 
be conceived to be, might exist in a spherical shape with a spherically symmetrical 
distribution of density. In the absence of gravitation there could be no question of 
instability. The effect of any local condensation would be to set up vibrations, and 
the frequency of the vibration of any spherical harmonic type would depend upon the 
elasticity of the material. If the resistance of the material to compression is suffi- 
ciently high the stability persists in spite of gravitation. There are thus two 
competing agencies : gravitation, tending to instability, and the elasticity of the 
material, tending to stability. In a general way it is clear that, as the elasticity 
diminishes, the frequency of vibration of any type also diminishes; and, if the 
frequency can vanish for sufficiently small elasticity, the planetary body possessing 
such elasticity cannot continue to exist in the spherically symmetrical configuration. 
The problem is to determine the conditions as regards elasticity in which the 
instability occurs. 

A grave difficulty presents itself at the outset. In the equilibrium configuration 
the gravitating planet is in a state of stress ; and, in a body of such dimensions as the 
Earth, this stress is so great that the total stress existing in the body when it vibrates 
cannot be calculated by the ordinary methods of the theory of elasticity. In that 
theory it is ordinarily assumed that the body under investigation is in a state so little 
removed from one of zero stress that the strain, measured from this state as a zero oi 
reckoning, is proportional to the stress existing at any instant. In order that this 
assumption may be valid, it is necessary that the strain which is calculated by means 
of it should be so small that its square may be neglected. Now if we apply the 
equations of the ordinary theory to the problem of a solid sphere strained by its own 
gravitation, and if we take the sphere to be of the same size and mass as the Earth, 
and the material of which it is composed to possess moduluses of elasticity as great as 
those of ordinary steel, we find that the strains may be as great as -J-, and thus the 
strains are much too great for the assumption to be valid. The initial stress existing 

* " The Stability of a Spherical Nebula," London, 'Phil. Trans. Koy. Soc.,' A, vol. 199 (1902), p. 1. 
t J. H. Jeans, " On the Vibrations and Stability of a Gravitating Planet," London, * Phil, Trans, Eoy. 
Soc.,' A, vol. 201 (1903), p. 1§7, Quoted below as " Jeans (1903)," 
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in the gravitating planet, the stress by which the self-attraction of the body is 
equilibrated, is much too great to permit of the application of the ordinary theory. 
The same difficulty presents itself in every problem concerning the elasticity of a 
gravitating planet, for example, in the problem of tidal deformation or of the stress 
produced in the interior by the weight of continents. In these problems the difficulty 
was turned by Lord Kelvin # and Sir G. H. Darwin! by taking the modulus of 
compression to be much greater than that of any known material, in other words, by 
taking the material to be incompressible. Their object was to determine the degree 
of rigidity which must be assigned to the Earth, and for that object it is permissible 
to turn the difficulty in this way. When the problem is that of gravitational 
instability this artifice cannot be adopted, because the whole question is that of the 
degree of compressibility which is admissible if the gravitating planet is to be stable 
in a spherically symmetrical configuration. The artifice adopted by Jeans (1903) 
consisted in annulling the initial stress by introducing an imagined external field of 
force to equilibrate the self-attraction of the planet. 

The problem thus posed is an artificial one, which may, nevertheless, throw light on 
the actual problem. When the initial configuration is taken to be one of uniform 
density, the analysis of the problem is of the same kind as that which presents itseli 
in the problem of the vibrations of an elastic sphere, a problem which has been worked 
out very completely by H. Lamb.J The determination of the effect produced by 
gravitation in lowering the frequencies of the various modes of vibration is reduced to 
a question of troublesome analytical computation. Jeans worked out the problem on 
the basis of the ordinary theory of elasticity, using the elastic constants X and /x of 
Lame. The constant //, is the modulus of rigidity, and the constant X is such that 
X + f/x is the modulus of compression. In the case of the Earth the values of these 
constants can be inferred from the observed rates "of propagation of the various types 
of disturbance which are perceived as earthquake shocks. He concluded that, when 
the proper values are attributed to these constants, the Earth must be held to be in 
a state far removed from one of gravitational instability ; but he suggested that, if 
the resistance to compression was at one time considerably smaller than it is now, the 
spherically symmetrical configuration would then have been unstable ; and he held 
that there are traces of the instability in the distribution of land and water on the 
surface of the globe. 

The actual problem differs from this artificial problem in the mode of balancing of 
the internal gravitation. Lord Rayleigh§ has proposed a method of meeting the 
difficulty as to initial stress. He proposed to consider the stress in the vibrating 

* See, in particular, Kelvin and Tait's 'Natural Philosophy/ Part II, §§ 833-846, Cambridge, 1883. 
t "On the Stresses caused in the interior of the Earth by the weight of Continents and Mountains/ 5 
London, 'Phil. Trans. Eoy. Soc./ 173, 1882, p. 187. 

J "On the Vibrations of an Elastic Sphere," London, 'Proc. Math. Soc.,' 13, 1882, p. 189. 

§ "On the Dilatatipnal Stability of the Earth," London, 'Proc, Boy. Soc./ A, 77, 1906, p. &86 % 
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gravitating sphere as compounded of two stress-systems : a hydrostatic pressure by 
which gravitation would be balanced if the sphere were in equilibrium, and an 
additional stress. He proposed to measure the strain, not from the unattainable 
state of zero stress, but from the equilibrium state ; and he proposed to take the 
additional stress to be determined in terms of the strain by those equations which are 
commonly used in the theory of elasticity. To simplify the problem he proposed to 
take the material in the equilibrium state to be homogeneous and the elasticity to be 
isotropic, so that the equations connecting the additional stress and the strain are of 
the same form as the ordinary stress-strain relations of isotropic elasticity. In 
justification of the proposed procedure he brought forward theoretical considerations 
founded upon the general theory of energy, and other evidence drawn from an 
interpretation of the experimental results in regard to the behaviour of elastic solid 
bodies. It is not too much to say that all the evidence there is, is just as strong in 
favour of Lord Rayleigh's proposed method as it is in favour of Hooke's law, in the 
sense in which that law is applied in the ordinary theory. The only objection which 
can be raised against the method, an objection mentioned by Lord Rayleigh himself, 
is that the body to be treated is certainly not homogeneous, and possibly not isotropic. 
When the proposed method is adopted, the density and the moduluses of elasticity 
must be taken to have their mean values. The justification for treating the values of 
these quantities at any point as equal to the mean values, is that it is advisable in 
the first instance to work out the simplest case. # 

In the first part of this paper the mathematical problem proposed by Lord Rayleigh 
is worked out ; and the conclusion is drawn that the effective moduluses of elasticity 
of the Earth, in its present state, are sufficiently great for a homogeneous spherical 
configuration to be thoroughly stable. The second part of the paper is devoted to 
developing the consequences of supposing that the elasticity of the material of the 
Earth was once much less than it is at present. 



Statement of the Mathematical Problem. 

2. We have before us a perfectly definite mathematical problem, which may be 
stated as follows : — A sphere of radius a, and of uniform density p , is in equilibrium 
under its own gravitation, and the stress within it is hydrostatic pressure of amount 
p at a distance r from the centre. When any small disturbance takes place, so that 

* It may be observed that the method advocated by Lord Eayleigh is the same, except for a slight 
modification, as that which was used in the second edition of my "Treatise on the Mathematical Theory 
of Elasticity, " Cambridge, 1906, in the discussion of the statical problem of a gravitating sphere held 
strained by external disturbing forces. The modification consists in the assumption, which was there 
made, that the material might be treated as incompressible. If this assumption is not made, the analysis 
becomes much more difficult. An earlier indication of the method will be found in a paper by J. Larmor 
Cambridge, <Proe. Phil, Soc.,' 9, 1898, p. 183, 
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the particle which was initially at (a;, y, z) is displaced to (x+u, y + v, z + w), the stress 
is specified by six stress-components X. x , Y y , Z z , Y z , Z x , X y , and these are connected 
with the initial pressure p„ and the displacement (u, v, iv) by the formulae 



X, 



•^ +x ^ + fi + t) +2 ^ 



3^ 

dx ■ 



-\ 



T.--ft+X(^ + ^ + ^V«C 3 ' 



dx dy dz 



dy' 



r? , x fdu dv , dw\ , dlV 



> 



Y.= 



(div dv 



fdu , Qw\ 



Y — /^ v a. ?!f N 



• • (1) 



J 



where X and /x are constants. It is required to form the equations of vibration, and 
to solve them, so as to determine the character of the modes of vibration and the 
equation for the frequencies, and, in particular, to ascertain the relations which must 
hold among the quantities X, /x, p , a in order that any frequency may be reduced to 
zero. We proceed to express this problem in terms of a system of differential 
equations which hold at all points of the body, and a system of special conditions 
which hold at all points of the undisturbed surface. 

3. In the equilibrium state the potential V at any point is given by the equation 



V = f irypo (3a 2 -r 2 ), 



(2) 



where y is the constant of gravitation. The equation of equilibrium is 

I r\ «— <> Vy » , « . » . a 



p dr dr 



or 



3j>0 4 



dr 



f rrypo r. 



. . . (3) 



(4) 



Since p = at the surface r = a, the value of p at any point is given by the 

equation 

p ^%7ry P() 2 (a 2 ~r 2 ). . ........ (5) 

When the sphere vibrates, the equations of motion are three equations of the type 



d 2 u 



8V 3X^ 3X^ 3Z* 



p a/, 2 p ax ' dx ' dy 



dz 



(6) 



where p is the density, and V the potential, in the disturbed state. In the left-hand 
members of these equations we may ignore the distinction between p and p . In the 
right-hand members we may put 

p = p (l— A), (7) 
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where A is the dilatation expressed by the formula 



A 



Further, we may put 



du dv J dio 

_— _ ~p _ , j- _____ e 

ax ay az 



» * • • a o • • » tOI 



v = v +w, . ...... 



(9) 



where W is the additional potential due to concentration of density at internal points, 
and to displacement of mass across the initial bounding surface. We may neglect 
terms of the type p A3W/3#. When we substitute for X^... from equations (1), and 
make these simplifications, the equation (6) becomes 



Po 



dt 2 



ox ax 



~ + (X + a) — + aV 2 U. 

dx dx 



On omitting the terms which cancel each other in virtue of equation (3), we have the 
first of the three equations (10) below. The remaining two of these equations 
are obtained in the same way. Thus we have the equations of vibratory motion in 
the forms 
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In addition to these equations we have the equation connecting the potential with 

the density in the form 

V 2 W = 4z7ryp A t . (11) 

The system of equations (10) and (11) are the differential equations of the problem.* 

4. Besides satisfying the differential equations (10) and (11), the additional 

potential W and the components of displacement u 9 v, w must also satisfy certain 

conditions at the surface r = a. Let U denote the radial component of displacement, 

so that 

JJr = xu + yv + ziv, . (12) 

and let U fl denote the value ofUatf = a. The potential W is that due to a volume 
distribution of density — p A, together with that due to a superficial distribution p Q JJ a 

* In the problem as formulated by Jeans, when the self-attraction of the body is balanced by an 
external field of force, the equations of vibratory motion differ from those which are obtained here by the 
omission of the terms such as ^7ryp i2 xA. In Lord Kayleigh's paper already cited, the equations given by 
Jeans are discussed in accordance with the analysis which was developed by Lamb in the paper on the 
vibrations of a sphere. 
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on the surface r = a. By the method of spherical harmonics we can, when W is 
known, write down the expression for the function W (0) which is the potential at 
external points of the same distribution. The surface characteristic equation gives 

'3W (0) \ /8W\ 
-5 — - -o-j = -47ry/)oU«, (13) 

This is one of the conditions which must be satisfied at the surface r = a. To 
obtain the other conditions which must be satisfied at this surface, we observe that 
the disturbed surface r = a + U a is free from traction. If l 9 m, n denote the direction 
cosines of the outward drawn normal to this surface we have three equations of the 

type 

ljL x -\-mJLy-\-nL x = 0, 

which hold at the surface r = a + U a . If in this equation we substitute for X^, ... 
from equations (1), we see that in the terms containing w, ... we may replace £, ... by 
the approximate values 05/7% y/r 9 z/r. The only term which does not contain u 9 , . . is 
the term — Zp arising from — IX. X , Now _p vanishes at r = a 9 and therefore at 
r = a + U a we have 

\0r/r = a 

to the first order in u, v 9 w. Hence in this term also we may replace I by xfr. On 
substituting for p from (5) we find that the equation 



r?ry/0o a?U+ - AA + 2//,— +^ /A _. + _) + - /jl _+ • = 
. r\ ox J r \ox oyj r \oz ox J 



must hold at the surface r = a. By an easy transformation this equation becomes 
the first of the three equations written in (14) below. The remaining two of these 
equations are obtained in the same way. The equations which must hold at the 
surface r = a are therefore equation (13) and the equations 



ju, ox or ft 
^ + £ (Ur)+ ,£_. +lT -,Ur-0.> (u) 

[x oz or ft 



These equations can be interpreted in the statement that the traction on the mean 
sphere is a pressure equal to the weight per unit of area of the material heaped up to 
form the inequality U fl . 

VOL. CCVII. — A. 2 A 
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5. We shall now suppose that the system executes a normal, or principal, vibration 
of frequency f>/2?r, or, in other words, that every component of displacement is 
proportional to the same simple harmonic function of pt. The equations of vibratory 
motion become three equations of the type 

( x + ^a^ + /xV2 ^ + p ° p2u + svyp*®* +/) °~^r = ' * * - • ^ 

where W satisfies the equation (11). The solutions of this system of equations (11) 
and (15) must be adjusted to satisfy the conditions (13) and (14) at r = a. These 
conditions can be satisfied only if p has one or other of a certain infinite set of values, 
which are the roots of the frequency equation. The problem of gravitational 
instability is solved when we find the conditions that one of the values of p may be 



zero. 



Solution of the Differential Equations by Means of Spherical Harmonics. 
6. We introduce the notation 

j£&- = h 2 9 ^ = * 8 , fa^!-6' 2 , W = 47ryp E. . . (16) 

A + 2/A [JL A+2/X 

The equations of motion (15) become three equations of the type 

(v . + i> + (*-l)f| + j^(«A + »§§)_0, (17) 

and the equation v 2 W = 47ryp A becomes 

y JL-i — ~ i—i. , . . « . » . • . * • I A O / 

in these equations A stands for 

du dv div 

__ — _ —4— »-— — ™ ■—♦— — 

ox cy oz 

By differentiating the left-hand members of the equations of type (17) with respect 
to 'x, y } £, respectively, adding the results, and simplifying by means of (18), we 
obtain the equation 

(v 2 + A 2 )A + 6s 2 A + s¥^- 0,. . . (19) 

or 

The method of solution of the problem is this :•— We seek first a solution of the 
equation (19) in which A has the form 

where (o n is a spherical solid harmonic of positive integral degree n, and f n is a 
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function of r which is such that r n f n is finite at all points within r = a> including the 
origin r = 0. We seek next to determine E in the form 

E = % [E w (r)(o w +F w ], (21) 

where E w (r) is a certain function of r which is such that r w E n is finite at all points 
within r = a, including the origin r = 0, and F n is a spherical solid harmonic of 
degree n. The equations of motion of type (15) then become three equations of the 
type 



ik 2 \ f 3 IF 

(v 2 +F)w + (^-i)S|^(/ w .^)J + ^s 2 ^S(/,.^) 



+ 'Ss 2 ^ 



h 



2 



dx 



in which we must have 



(E n .a) n )\ + j2-s 2 t 



/3F. 



w 



\3x 



= 0, . 



du dv dw ^ , r \ 



3x 3?/ 82 



It appears on trial that we can obtain a solution in which 



where u u v u v,\ satisfy the system of equations 

(V+A")«i + (^-l)s{^(/..«-) 



£2 



^ 



A s 



+ ^s 2 $ 
h 2 



dx 



(E„ .<o n )\ = 0, 



V 



• J 



3^1 . 3^i , 3^x ^ / /. x 
also m 2 , v 3 , w 2 satisfy the system of equations 



J 



% 
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8x 



v/ 2 • • » 
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3^2 . 3^2 , 3^2 a . 

~r ^ r -^ — — U ? 



8x 3y 82; 



y 



and t%, v 3 , ?# 3 are a complementary solution of the system of equations 

(v 2 +*>3 = 0, (v 2 +F)^ 3 = 5 (v 2 +F)t^ 3 = 0, 

3w 3 



3^3 + 3^3 + 

8x dy dz 



> 



0. 



. (22) 



(23) 



u = Ux+Ug+Ug, v = Vx + v 2 + v 3 , w = Wi + w 2 +w 3 , .... (24) 



. . (25) 



(26) 



• • 



• W 



7, The sets of functions u u v u %v x and u 2 , v 2 , w 2 can be any particular solutions 

2 A 2 
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of the systems of equations (25) and (26). It appears on trial that u u v u iv x can 
have the forms 



u x = t 



)0) 



v, 



$ 



Uh 



J?.(r)^ + Q.(r)t 



cz 



2n + 3 5_ 

dx 


( 0) n 

\ r 2n + l 


2n + 3 U 1 




,2w + 3 ^ i 


/ 0) M 



v 



3z \r 



2» + l 



. . . . (28) 



where P„(r) and Q„(r) are certain functions of r. Also it is clear that u 2 , v 2 , w 2 can 
have the forms 



3s 2 ^ /3F 



n 



v> ~ - W i ^ 



w, = 



3s 2 ^ 3F„ 



A 



3z 



(29) 



Further, the forms of ^ 3 , v& w 3 are known from the analysis of the problem of the 
vibrating sphere which is free from gravitation. We have 



u 3 = t 
where 



>< fc )(»|?-*|?) + *-( h -)^ - „Il*-< tr >' A "'4(^*.)} ' (30) 



*" (fo ' ) = *&^s)^& i) • (31) 



^ n and <j) n are spherical solid harmonics of degree n> and the expressions for v s , io 3 are 
to be obtained from the expression for u 3 by cyclical interchanges of the letters x, y, z. 
It appears that to a single term f n co n in the expression for A there corresponds a 
definite term ¥ n in the expression for E. Further, when we form the boundary 
conditions, it appears that the terms of it 3 , v 3 , iv 3 which contain ^ represent a free 
vibration, and the frequency of this vibration is determined by the same equation as 
if the sphere were free from gravitation. Also it appears that to any term f n o) n in 
the expression for A there corresponds a definite function <f> n in the expression for 
u 3 , v 3 . w 3 . The solution expressed by a single term f n a> n of A and the corresponding 
terms of (u u v u iVij, (u 2 , v 2 , w 2 ) and (u 3 , v 3 , w 3 ) determines a normal mode of vibration. 
We shall therefore omit ^ and reduce all the summations to typical terms. 

8. If in the equation 

3A 



(V 2 + A 2 )A-f6s 2 A + s 2 



dr 



(19 bis) 



we put f n (r) . o) n for A, we find that f n (r) must satisfy the equation 

dr r dr r dr \ dr / 



or 



n 



dr 



,2 



+ \ 



2 ( n + 1 ) +s 2 r\^ + {h 2 +(n+6)s 2 }f n = (32) 

v J dv 
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This equation is a linear differential equation of the second order ; and the forms of 
the coefficients show that the point r = is a critical point, and that there is no 
other critical point at any finite value of r. If we seek a solution in series of 



the form 



fn = <W W + Ol* 



• m+1 + a 2 r m+2 +..., 



we find the " indicial equation " 

m (m— 1) +2 (n+l)m = 0, 

from which either m = 0, or m = — (2n+l). We must take the series for which 
m = 0, because r n f n must be finite at r = 0. Further, the form of the equation shows 
that this series contains even powers of r only. We assume, therefore, for f n 

the form 

f n = A[l + a 2 i* + ap* + . .. + a 2/f r 2K +...] ? 

where A is an arbitrary constant, and then we find the sequence equation 



a 



2k + 2 



(2K+2)(2K+l + 2n+2)+a 2K {h 2 +(n+6)s 2 +2KS 2 } = 0, 



or 



#2k + 2 — 



a, 



2k 



h 2 +(n-i-6 + 2K)s' 



Hence we have 



(2k+2)(2ic+2w + 3) 



fn — A 



1 _ h l+ ( n + Q ) s \ *\ IK+ 0l±6) s 2 }{h»+(n + 8)s*} 4 _ 
2.(2w+3) 2.4(2w + 3)(2« + 5) 

, y {^+(n + 6)^}{^+(n + 8)^}...{^+( W + 2 K +4)^} , K 1 , . 

V ■' 2.4...2«r(2» + 3)(2n+5)...(2»+2K+l) K } 



The series is convergent and represents the function^ for all finite values of r. 
9. We must next determine the function E n (r) from the equation 

V 2 (E B a>„) =f n (o n , 



or 



d 2 E n + 2(w+ l) dE„ = - 
a?r 2 r e/r ' 



(34) 



or 



^{,-f, + ( 2 » +1 )E,}=,/, 
We introduce an intermediate function 6 n (r) by the equation 



^ = r^» + (2n + l)E.. 



(35) 



Then 



#» 



-J* 



e£r 



^C+A 



1^2 



r ~_ 



A»+(n + 6)*» {A 2 +(n+6) g 2 }{fr 2 + ( n+8)s 2 } . 
2.4.(2n+3) 2.4.6(2w + 3)(2w + 5) 
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where is constant of integration. Then we have 



d 

dr 






and therefore 






c 



c 



r fc+1 + 2»+l 



+ A 



/v»2 



{h 2 +(n+6)s 2 } r 4 

2 . (2n+ 3) 2.4 (2^+3) (2^+ 5) 



+ ... 



where (7 is a constant of integration. Since r n E n is finite at r = 0, the constant C 

must be zero ; but the constant C is in our power, and we may choose it in any way 

that is convenient. The term contributed to E by C is (2n-f l)" 1 Cco n , which satisfies 

Laplace's equation, and therefore any change in the chosen value of C is equivalent 

to borrowing a term of F n to make up a term of E n o) n . 

Now the series 

1 _ h 2 +( n + 4) s 2 r2 {h 2 +(n + 4)s 2 }{h 2 +(n + 6)s 2 } ^_ 
2/(2n+i) r " 2.4:(2n+lj~(2n+3) ? *" 



satisfies the equation 



d " + ( 2n + # r \ * + {h* + (n + 4) s 2 } 



d? 



r 



dr 



1 



h 2 +(n + 4:)s 2 2 



2 



(2w+l) 



and therefore, if we take for C the value 



C = 



2n+l 



h?+(n+A)s 



A 

2 ' 



the function n satisfies the equation 



0; 



d% 
dr' 



'2n 



dd. 



+ 1 — + s 2 r)-^ + {h 2 +(n+4:)s 2 }d n = 0. 



. . (36) 



We shall choose this value for C, and thus we shall have 



0» = A 



2 w+l i r2 _ h 2 +(n + 6)s 2 

h 2 +(n + 4)s 2 2 2.4(2n + 3) 



r*+ ... 



and 

E. = A 



/ y +1 {h 2 +(n + 6)s 2 } {h 2 + (n+8) s 2 } ...{h 2 + (n+2 K +2) s 2 } , 3k ' 
1 ' 2A...2 K (2n+3)(2n+5)...(2n + 2K-l) r "*'. 



,.2 



^ + 



{h 2 +(n+Q) s 2 } r 4 



_ h 2 +(n + A)s 2 2(2n + S) 2 . 4 . (2n + 3) (2w+5) 



I « ».._• 



The function E„ satisfies the equation 



dr 



n , , 2(w+l) | g2 1 cCE, 
r J c£r 



^ ' 2.4...2«r(2n+3)(2n+5)...(2n + 2K+l) r "**. 






(37) 



(38) 



+ {h 2 +(n+4)s 2 }E a = (39) 
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It will be convenient presently to have observed that the equation derived from 
this one by differentiating the left-hand member with respect to r can be written 



d' 



dr 



- + i 



2_(«±2) +sV l d + „ s+ , ^y, 

r j dr 



1 dE n 



r dr 



0. 



. (40) 



10. The forms of u 2 , v 2 , w 2 and u 3> v s , w% have been put down and it remains to 
determine u u v u n\. We have a system of three equations of the type 



(V 2 +F) Ui + ('£- 1) {f*g + 1 & x<o n ) + p%x<o 



Iv 



dx r 



+ 



h 2 



i w, 



IdE. 



i' 2 (E.,^-" + ;^-"xft. 



•n 



dx r dr 



n 



0. . 



• • 



(41) 



We express xco n in the form 



xo) 



ry* A 



>ft) 



n ™2n + 



ft) 



n 



2n+l [ dx 



3^ \r^ + \ 



(42) 



and then the above equation becomes 



(V 2 +FK = - 



Oft) 



tt 



3# 



A 2 



+ 



1 



^V n 



n \ i 2 

2^+1 ofr / A 2 " 



,™« 



,2rc,+ 



-/ n +3E n + 



3r c£E n v 



2n+l dr 



,,2n + 3 



ft) 



ra 



2»+l 



dx 



2n+l 



,r 



h 2 






.(43) 



We seek solutions of the system of three equations of this type in the forms of 
the type 

(44) 



v- 



P.^+Q,*** 8 



3a? 



dx \r 2n+1 ' 



in which P„ and Q„ are functions of r. We find 



(V 2 +F)«, 



n 4- —. n -i_ PP 



d/r 2 



<fr 



ra 



Oft) 



ra 



8^ 



+ 



' d 2 Q n + 2(n + 2) dQ % 
dr 2 r dr 



+ k 2 Qn 



,2?i + 3 



ft) 



w 



3^ \r : 



2n+l 



Hence the assumed forms can be adjusted to satisfy the equations if V n and Q n 
satisfy the linear equations 



n 4- _ _f_g 4_ 7/< 2 p 



d!r 



r dr 



n 



d 2 Q n + 2(n+2) dQ n + F ~ 



,'V 
h 2 ' 

A/ ^ 

h 2 \2n+ 



'¥ 



!)/-+; 



ra 



2n+l dr 



r 



- /, + 3E n + 



3r <iE n 



2n+l d!r 



• * 



(45) 



dr' 



r 



dr 



2n+l [\h 2 



1 )rdr + h' 8 V n+ r~dV)\- ^ 
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The right-hand members of these equations can be simplified by introducing the 
functions 6 n and r~ l dE n /dr. We have 



fn + 



r 



df n __ 1 d9 n r d (1 d0„ 






2n+l dr r dr 2n+ 1 dr \r dr 



2n + 1 \ dr 2 



d 2 6 n 2n d0, 



ii 



dr 



and 



/y»* 



2n + 



T /.+ 



3E W + 



3r dE n 



2n+t dr 2n+l \ d? 



r^L n + 30, 



n I ? 



and thus the equation for V n becomes 



n _i_ zUl „j_g X. PP 



dr" 



r dr 



n 



l 



'd\ 2ndd. 



n 



1 k 2 [d 2 e» . 2n dd 



d6» 



2n + l \dr* r dr J 2n+ 1 h 2 \d^ 



n + — ^ + 6 >V ^ + 3,S 2 tt 



r 



dr 



dr 



1 /cP0 n , 2n dd n 



- 



+ 



2^4- 1 \c£r 2 r dr 



-+W 



2w+l 



^-^ 2 {A 8 + ^ + 4)^^ + 3^^ 



by the equation (36). Hence we have 






n 



2n+\ 



4-P' 



(47) 



where P'„ satisfies the equation 



d 2 P' n , 2n d?' n up, _ w+1 P ™ 

_j_ — _ ~y- fh J. n — ~ ~ T7i « v «,« • » 



dr 2 



dr 



2n+l h 



2 « v n« 



(48) 



Again, we have 



I 
r dr 



n 



I A J^ 2 Ey 2(n+l)dE 
r dr 1 rfr 2 r dr 



?i 



i d 



r dr 



d [JldE n \\2(n+l)dK 



dr 



x dr 



+ 



dr 



and 



d 2 fldE H \ 2(^+2) d /ldE n 



dr 2 \ r dr 



r dr \r dr 



. 3 dE„ d 2 E, , 2 (n + 1) dE n , 3 dE n 



r d 



dr 



r 



. -4— -«, 

dr r dr 



= r A A ^ + 2 ( n + 3 ) ^ . 

dr \r dr J r dr 
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and thus the equation for Q„ becomes 

d 2 Q n 2(w+2) dQ n pn 
ar r dr 



2n+l 



cZ 2 /l dE n \ 2(n + 2) d (I 



n 



dr 2 \r . dr 



r dr\r dr 



+ 



2n+l¥ 



(P_ (I dE n \ 2 (n+2) d^fldE 
dr J 



dr~ \r 



r dr \r dr 



2 d (\ dE n x 
+ s 2 r l 



dr \r dr , 

+ 2(n+3)s 2 (±zp. 

\r dr 



+ 



2n+l 
2^+1 



d 2 (I dE n \ 2(^+2) d /l dE B \ , 2 /l dE 



dr* \r dr / 



n 



r dr \r dr 



r dr 



,r 



F ( 1 ^pA _ *! {#■+ (n + 6) s 2 } fi —^ + S 2 (w + 3) s 2 f- ~ 
T rf?" I ii \r dr h \r dr 



n 



by the equation (40). Hence we have 



Q. = -^-^=+Q' 



2n+l r dr 



• • ■ 



. . . (49) 



where Q' n satisfies the equation 



d 2 Q! n + l{n±2) dQn +k2 Q, _ n Jc 2 o2 1 dE 



dr' 



r 



dr 



-o s' 



2n + 1 h 2 r dr 



. . . (50) 



11. To find the form of P ; w , we assume 



F * = £jT[ h 2 S * A ^° +p ** +P ^ + ~^ 



• • • 



(51) 



Then equations (48) and (37) give 



Fp +2(2n+l)p 2 = - 



2n+I 



k 2 p 2 +i(2n + 3)p i 
Fpi+6(2n + 5)p 6 



k 2 p 6 + 8 (2n + 7) p B 



h 2 +(n + ±)s 2 ' 



x 

2 > 



_ fr 2 +(w+6)s 2 
2.4.(2n+3)' 

{A 2 +(n+6)g 2 }{F+Qi + 8)s 2 } 
2.4.6(2» + 3)(2n+5) 



k% K + (2 K +2)(2n+2K+l)p 2K +. 



/_!).+! {A 2 + (n+6)s 2 } {/a 2 + (^+8)s 2 }...{4 2 + (ra+2/c+2)s 2 | 
' 2.4...2ic(2w+3)(2n+5)...(2rt+2K-l) 
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By these equations p is left arbitrary, and p 2 , _p 4 ... are determined when p is 
chosen. As we need a particular integral only of the equation for P'„, we may choose 
p in any way that may he convenient. We shall put p a = 0. Then 



P* = - 7 



1 



21 ' 



2{h 2 +(n + 4:)s 2 } 



P4 = 



+ 



t 



2.4(2n + 3) 2.4.(2«, + 3)/i 2 +(» + 4)s 2 ' 



P«= - 



2.4.6(2n+3)(2w + 5) 



h 2 +(n + 6)s 2 + k 2 + 



¥ 



h 2 + (n+4)s i 



Ps 



2.4.6.8(2n + 3)(2n + 5)(2n+7) 



{h 2 +(n+6)s 2 }{h 2 + (n+8)s 2 } 
+ F{h 2 +(n+6)s 2 } + Jc i + u 



k 



.6 



h 2 +(n + 4)s i 



To find the form of Q' n , we assume 



W> n 



n h 2 



2n + 



Then equations (50) and (38) give 



k 2 q + 2(2n + 5)q* 



1 



2n + 3 



7 9 * /« h\ h 2 + (n + 6)s 2 

k 2 q 2 + 4 (2n + 7) q± = - —,~ — - oWo 7 , .x , 
* v /:2 2.(2n+3)(2n + 5) 



Pg 4 + 6(2n + 9)2e 



{h 2 + (n + 6)s 2 }\h 2 + (n+8)s 2 } 
2.4(2n+3)(2n + 5)(2?i + 7) ' 



> 



> 



j 



(52) 



-r 2 s 3 A(g +f/+^* 4 +...) (53) 



72 /o ovo ,o L «\' / v {^+(»+6K}{/ t 2 + (n + 8) S 2 }...{^+(n+2 K +4) S 2 } 

g gfc +(2«+2)(2n+2,e+ %»+,= (-)■ •! ^^^j (2n + 3) (2n + 5)...(2n + 2ic+3) ' 



As before, g can be chosen arbitrarily, and then q 2> q i9 ... are known. 

that if we put 

_ 2p 2 _ 4p 4 ^ y _ ( 2k + 2)j3 2k+2 



We observe 



ffo 



2n+3 



g 2 = 



277,+ 5 



?* = 



2^ + 2k + 3 



the sequence equations for the g ? s are transformed into the sequence equations for the 
p's, beginning with the equation containing p 29 jp 4 . We shall therefore choose q to 
be (2n + S)" 1 2p 2> and then 



Sy n 



n ¥ 



\n+l H 



s 2 A 



2p 2 , 4p 



+ 



4 n,2 



?■* + 



6p 



6 n „4 



2n+3 2n + 5 2n+7 



r +... 



? * * 



• (54) 
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This choice of the q's amounts to subjecting the functions P' B and Q f n to the 

equation 

clV 



dr n+1 



dr 



(55) 



To see that this equation is compatible with the differential equations (48) and (50) 
for P' n and Q' m we observe that 

du n __ 1 d f 2n +sf _1 ^-E^ 



rf? 



J2n + 1 



and that from the equation (48) we can form the equation 



d 2 2(n+l) d ^ 72 2(n+l)l 
dr r dr r } 



r 



2n + 



,dP 



/ \ 

n 



n+1 h 2 o d 



dr 



2n+lh 2 dr 



r 



2w + 3 



IdEn 



x dr 



while from (50) we can form the equation 



dP _ 2 (n+1) d_ p 2 (n+1) 



dr 



r 



afr 



r' 



> < 



d_ 

dr 



(r 2 » +3 Q' re ) 



n P^^J^,/! dE 



w 



2n + 1 A 2 ctfr 



T (ir 



12. We have still to satisfy the condition 



ox oy oz 



(56) 



If we form the expression in the left-hand member from the expressions of the 
type (44) for u u v l9 w l9 we have 



du x dvx dtVx 

— . _j_ — _j_ 

dx ' dy dz 



\r dr ™ 2n+2 ~ 7 ~^ ' ' 



r 



2n + 2 dr 



By means of the formulae (47) and (49) for P w and Q n , the coefficient of co n in the 
right-hand member of this equation is transformed into 



2n+l 



~n d6 n n+1 d_ f ^ 2n+2 



r dr 



r 



2n + 2 



dr 



dr /_ 



+ 



n + 1 

„2n + 2 



n 



,2n + 



1 \Ajjl n Cv 



n+1 



- it (r 2 " +3 Q'«) 



dr dr 



The first term is equal to 



—-—[nf n + (n+l)f n ] or f H9 



and the second term vanishes identically in virtue of equation (55). It follows that, 
with our choice of p and q 0i the equation (56) is satisfied identically. 

13. We have now completed the determination of the forms of u 9 v, iv in terms of 
the spherical solid harmonics a> m F tt , <j> m and of certain functions of r, viz. : f n9 n9 E n , 

2 B 2 
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P'n> Q'n> $n> Various relations between these functions have been noted incidentally. 
It will be convenient hereafter to have noted the following properties of \jj n (Jcr) : — 



k 2 r 






¥r 



\b (hr) = LZli! i 1 — - 

^ v ; 1.3.5...(2n+l) I 2(2w+3) ' 2.4(2n+3)(2n+5) 



(57) 



d* 



<i (/cr)' 



+ ^ ti) rfffi) +i }*-( fe ')= o - < 58 > 



£r^|^ , . (59) 



Adjustment of the Harmonies to Satisfy the Boundary Conditions. 

14. In order to express ¥ n in terms of o) n and <f> n we use the condition that 
47ryp (E n a) n + 'F n ) is the potential at points within the sphere r = a of a distribution 
of volume density within the sphere and of surface density on the sphere. The 
corresponding external potential is 

^7Typ (a/r)*** 1 {E n (a) co w + F„}, 

where E n (a) is the value of E w at r = a. The surface density is p U a , where U a is 
the value at r = a of the radial displacement U. Hence we have the equation 



dr 



ax 

L.\? / 



2^ + 1 



{E n (a)co tt + F n } 



— — (E n co n + F n ) = - U. 



which holds at the surface r = a ; it gives 



2n+l 



cZE, 



2n + l 



r 



E n +^-")co n +^^^F 



dr 1 



r 



w 



U. 



«• 



r — a 



N 



ow 



n 



3r*s 



n 



U = l t P l ,_(„ + i)rQ ll U,- ^gF K + ^{^(k^+Fr^ilcr)}^ 



It follows, that the equation 



2n+l , 3ns* 2 

r 



a 



F = CO 

a/^ 2 J " B 



r 



P w -(^+l)rQ w 



dE n __ 2n+l 
dr r 



E n 



r = a 



n 






holds at the surface r = a. Since this equation connects the values at r = a of three 
spherical solid harmonics of the same degree n, it holds for all values of r, and gives 
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a generally valid expression of ¥ n in terms of a> n and $ n . By means of equations (35) 
and (59) the equation becomes 



im 



h 



n 



3ns 2 + (2tt+l)ft 2 ^*" nF|>+ i n+1 ) a2 ^}^- ( °n+n(2n+l)xjj n (ka)<l> n l . (60) 



15. The three remaining conditions which hold at the surface r = a are expressed 
by equations of the type 



yt, — 2 )xA+-z-(rTJ)+r- w+^^rU = 0. . . . . 

h 2 J dx K ' dr nr 



• (61) 



Every term of the left-hand member can be expressed in terms of the spherical 

solid harmonics 

dcjy 9 



uo) n r 2w+3 v I Wn 



._— „ l | — . ,,y, ,_ 

We have 



dx ' 



tX/L* — — JUG) % • / 



8a? \r 



2w+l /» 



3a; * 



n r 2?i + 3 & { Y 



2n+l 






m2 



2n + 



1^ 



JO) 



„2w+3 



?z 



3x \r 



f-^l ^n 



dx 2n+l dx\7 



2n+l 



Also 



rV = {nP w -(^+l)r 2 Q w }^-7i(2^+l)^(^)(^ 

3 ns 2 



n 



m j m* 



and therefore 



_— ™^— jyp [{0 n -nP n + (n + l)a 2 Q K } r = „.o> w + n(2n + l)^(te)^] ; 



3x 



(rU) = 



Again, 



3^ 
r ^ 

or 



nP w -(^+l)r 2 Q re 4- 



w(2ft+l)i/f ft (&r) 



3ns s 



3ns 2 +(2n+l) 

3m* 2 



^~ 2 {0 a -wP R + (n+ l)a a Q»} r 



« 



)<w 



w 



3x 



3?is 2 +(2^+l)/ 
dm 



JT-p^(te) 



3^ 



+ ^ { «p,-(» + i)^,}^ 



dx 



n ^2»+l J^( _ W w 



-4/yj \ /y»<"W« A ^ 



d\k n (kr) 



dr 



dx l dx\r 2n+1 



tyt 



d? 



dr 



^ + (n-2)P„ 



i<y 



n 



dx 






d 



3 / 



r fr {7 * n+SQn) - {n+Zyn+3Qn ldx\ 



0) 



n 



„2n+l 



+ 



3(n-2)s' 



Sns 2 +(2n+l)h 2 



!0 n -nY n +{n+l)a 2 Q n } r ^+n{2n+l)x}; n (ka) d< f > 



n 



I 



dx 



dx 



+ w d^M + {n -2)^_ 1 (kr)} d 4 



dx 



n 



n+l 



M 



d 



r ±. {7*+ B $ n+1 (kr)} - (n + 3) r 2n " ? % +l (kr) 



d J j>n 

/dx\r* +1 J' 
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xrXJ 



nP„-(w+l)r 2 Q M 

3ns 2 






3ns 2 +(2n + l)lr 



{0 n —nP n + (n+l)a 2 Q n } r=: 



rp' 



2 



-n(2n+l)\ifi n (kr) 



3ns 



.2 



3ns 2 +(2n+l)ti 



$ n {ka) 



/y»* 



2n+l 



_2n+l 

dcf) 

dx 



0)91 -^+1 



O) 



\ n 



n 



dx 



n _J2n + l 



dxyr^ 
dx\r 2n+l J\ 



»+i 



16. The equations which hold at the surface r = a can be arranged in such 
forms as 



A " 4-B r 2n+z ° I *°n 



dx 



dx\7 



,,2w+ 



1 ) +C W 






3 ^(^)^ 0? * # • (62) 



in which A ft , B w , C n , D n are certain functions of a, viz. : 



•AXjt 






j/n + 



3ns s 



0n+ l J? n Z l)h \,2 {^n- (»+l)«U} 



+ 



+ 



+ 



« c£ 



3«s 2 +(2n+l)A 2 " 3ns 2 +(2n+l)h 

dP. 



2n+ 



--r- {n¥ n -(n+l)a 2 Q n } +a-j^ + (n-2)P» 



- — p-^r- ^ W9 [0„ - ftP w + (ft + 1 ) a 2 QJ 
3ns 2 +(2n+l)h 2L v 7 ^ J 



s 2 /c 2 a 2 



(2w+l)# 



3ns' 



*• + o l^/o^XM {"P." ( W + 1)«U} 



3ns 2 +(2w + l)P M 3ns 2 +(2n + l)h; 



, ■ (63) 






» 



¥ 2\-J— f. ^— - ~ {nP n - (n+ l)a 2 Q„} +a^ +nQ„ 

+ 1 ada x ' da 



.h 2 J2n + 
s 2 F 



1 /»" 



(2w + l)/i 2 



2m 
3%s 2 



0» + ; 



(2n+l)h l 



_3ns 2 +(2»+l)/i 2 n 3ns 2 +(2w+l)/i, 



{nP s -(n + l)a 2 Q„} 



, . (64) 



c 



n 



—n 



(2n+l)xp n (ka)+a- r xfj n (Jca)-- ~ ; o ^ti\^ ^( fo ) 



c£a 






3ns 2 +(2n+l)h 

W+(2n+l)y ^ (to) + °55*- 1 (* a ) + ( w - 2 )*-x(* a ) 

_ w(2n + l)s 2 &V , /t \ 

3w 2 +(2«+lU 2<M '' 



• * * • • 



(65) 



D 



n d 



y\f n (ka) 



n 



d 



a da —-/ n+ l k \ a TJ^ k ^ + ™M te ) 
n(2?i+l)s 2 A: 2 



+ 



37i5 2 +(2^+l)^ 



$ n (ka). . 



• • 



»•*••«• 



• (66) 



In these formulae a is supposed to be written instead of r in the expressions for 
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17. We may express A w and B„ in more convenient forms by the use of the 

identities 

7iP n -(n+l)a 2 Q n = n -(n + l)l& n +n¥' n -(7i + l)a 2 Q' n , 



da 



a Jn> 



d 2 E n | 2(n+l)dE H = , 
da 2 a da 



a d ^ = e n -(2n + l)E n , 



a ^91» +(2n + 3) Q' B = -5L ± ^ . 
da n+1 a da 



n 



We find 



• = Pi£i« + "-> + » 7GT' +2 (»- 1 '(^T + P '" 



+ 6 3t-+ 1 ( ) ^T^ {( " +i)E -~ np, - +(n+i)gSy - > - • • (w) 

+ 3^ H- (2n + 1) A 2 {(n+ 1} E ' ~ nY ' n+ (W+ 1} a2Q/ " } - * ' ' (68) 

We may also express G n and D„ in simpler forms by the use of the equations 
connecting the \fj functions. We find 

C„ = -2(n-l)^[fc +1 (l») + z Jl n { ^ 1)hi |' £*.(*») + 3 ^ +( ^ +1) y *■(*») 

, (2w + lUW , /7 x ,_. 

^ { ^( W8B+ ^ lH+ _^^ Mfa) ], .... (70 ) 



D„ = # 



27ite Frequency Equation and the Condition of Gravitational Instability. 

18. Exactly as in the problem of the vibrating sphere which is free from gravitation, 
it follows from the equations of type (62) that we must have at once 

A„G> W 4-Cy>„ = 0, and B w a> w 4-D w ^ w = , ...... (71) 

and the frequency equation is of the form 

AJD w -B w C„ = (h ......... (72) 

The forms of all the functions which enter into the expressions of A w , B WJ O m D n have 
been determined. 

To investigate gravitational instability, we have to determine the conditions which 



192 PEOFESSOE A. E. H. LOVE ON THE 

must hold in order that the frequency equation may be satisfied by p 2 — 0. When p 2 
vanishes, h 2 and k 2 also vanish, but the quantity k 2 /h 2 , which is (X-f-2jn)//x, has 
a determinate value. We may not, however, obtain the result which we seek by 
first replacing k 2 /h 2 , wherever it occurs, by (X + 2/x,)//x, and then putting h 2 and k 2 
equal to zero wherever they occur, otherwise than in the ratio k 2 /h 2 . This pre- 
cautionary statement is necessary because it appears from the formulae (69) and (70) 
of § 17 that G n and D n both vanish if h 2 and k 2 vanish. Thus we ought to regard 
the equations (71) as being equivalent to the equations 

K<»n + (G n k~ 2 ) k 2 $ n = 0, B w Gi n + (D n k~ 2 ) #ty n = ; 

in other words, we ought to remove a factor k 2 from the equation AJD tt — B n C„ = 
before putting h 2 = and k 2 = 0. An exceptional case occurs when n = 1. In this 
case G n k~~ 2 vanishes when h 2 vanishes, and it will appear that A ?l also vanishes with h 2 , 
and the equations (71) ought to be regarded as equivalent to 

( AJi' 2 ) o, + {OJi- 2 k- 2 ) (k 2 ^) = 0, B^ + J) x k- 2 (k 2 ^) = 0, 

and we must remove a factor h 2 k 2 from the equation AxDj — B^ = before putting 
h 2 = and k 2 = 0. When we proceed in this way, the equation A tt D tt — B n C n = 0, 
with the appropriate factors removed, and with h 2 and k 2 put equal to zero after 
their removal, becomes an equation to determine s 2 a 2 , or ^Tryp 2 a 2 /(\+2fi). If the 
equation has a real root, the value so determined for s 2 a 2 gives a value of X 4- 2/x for 
which instability can occur. Since k 2 <j> n is a finite multiple of co n when X+2/x has any 
such value, it is certain that the homogeneous spherical configuration really is unstable 
for such values of X+2/x. If the value of X~h2/x belonging to the body is but little 
greater than the critical value, the equilibrium is practically unstable ; for a large 
displacement takes place if the sphere begins to vibrate according to the type 
specified by the degree n of the corresponding spherical harmonic function. For 
practical stability it is necessary that the value of X+2/a should be well above any 
critical value. The equation which yields the critical values contains the constant 
(X + 2fi)/[jL as well as s 2 a 2 . It will be convenient to write 



fl K A 



v =t^ = i (73) 



The value of v Cannot be negative, nor can it be greater than f . If the Poisson's 
ratio {X/2(X+/^)} of the material is positive, v cannot exceed J. If the modulus of 
rigidity /x were very small in comparison with the modulus of compression X-ff/x, 
v would be very small. If the velocity of propagation of waves of dilatation were 
twice that of waves of distortion, v would be \. This appears to be the most 
appropriate value to assume in the case of the Earth (see § 40 below). Since it is 
improbable that the ratio of the rigidity to the modulus of compression of the Earth 
has diminished since the date of consolidation, it will be sufficient for our purpose to 
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examine the two cases in which v = and v = |. We have now to discuss the 
conditions of gravitational instability in respect of the values 0, 1, 2, ... of the number 
n which specifies the type of vibration. 



Instability in Respect of Radial Displacements. 

19. The case in which n = is the case of a sphere vibrating radially. This case 
is not very easily included in the foregoing analysis, and it is very easy to investigate 
it independently. Let U denote the radial displacement. Then IT is a function of r, 

and we have 

^tt 2/tt 2-tt a dJJ , 2 IT 



r 



r 



dr 



r 



We go back to the equations (15) of the type 



(X + /i,)^- + /xV 2 u + p p 2 u + i7ryp 2 xA + p - — = 0, 



3# 



3^ 



where W, the additional potential, is a function of r. This equation is 



x d /dU . 2U\ 



r dr\dr r j 



i a? 



#+2.£V^) + 



^ r 2 r c l r \ r 



r ar\r j J r 

fdU , 2U\ , a; dW 



Now 



+ s 2 (\+2n)x ^ + — +Po-^ = 0. . (74) 

\ar r j r dr 



d 2 W x 2dW . fdU 2tT 



dr 2 r dr 



dr 



r 



and therefore we may write 



dW 
dr 



= 47ry/)oU + R 



where 



— - + -R = 0, or Rr 2 = const. 
dr r 



Since cEW/dr is finite at r = 0, we must have R = 0, and thus equation (74) becomes 
after division by (\+2/x,) xjr 



^ + 2^_^ +g ,/ ^ + 2U) + 3«"U+A 8 U = 0,. . . . (75) 



o!r 2 r dr 



where s 2 = f7ry/> 2 /(\+2fi,) and A 2 = p 2 p /(^ + 2ju.). This equation can be solved by 
means of a series, which is convergent for all finite values of r, in the form 



U = A 



VOL. CO VII. A. 



> _ h2 + 6s * r s , (ft a +68 8 )(ft»+8s 8 ) , 
.3 2.3.5 2.4.3.5.7 

/_ y (h 2 + 6s 2 )(h?+8s 2 )... {V+(2K+4 : ) s*} r2K+1 
K ' 2. 4. ..2k. 3. 5...(2k + 3) 

2 o 



• (^) 
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where A is an arbitrary constant. The second solution of the differential equation 
for U becomes infinite at r = 0, and thus the above is the most general form for U. 

20. The condition that the surface r = a + JJ a is free from traction can easily be 
shown to be the condition that 



4iryp VU+(X+2u)^ +2X— = 

dr r 



at r = a> Hence we have 



(f). + < 2 - 4 " +sW >v = ° < 77) 



where v = /x/(X + 2/x), so that 2X/(X + 2/x) = 2 — 4f. The frequency equation is 
therefore 

l(l + 2-4F + sW)^|~^a 2 (3 + 2-4^ + 5 V)+... 



Z . o . O 



( _ y .(^_^) J^+8^)...{^+(2K +iK} 2K( 22) = 

K ' 2.4...2k.3.5...(2 k +3) 1^+^ 4v+sa;... u. . 

The condition of gravitational instability is obtained by putting h 2 =0. It is 



(78) 



^(Z-iv + sW) 



Qs 2 a 



2„2 



^j . o . 



(5 — 4*> + sW)+... 



(-)' 6.8...(2*+4)s 2 "a 2 « _ ( 2lc +3-4 y + ^a a )... = 0. . (79) 



2 . 4. ..2k. 3 . 5...(2/c+3) 
21. The coefficient of s 2k q 2k in the left-hand member of (79) is 



8 1 3.5...(2/<+3) v ; 



or 



(-)' 



(2;<+3) 
_L3. 5...(2k-1) + 3. 5...(2fc+l) 



2k(2<c+2) 
3. 5.~.(2k+1)J' 



V < 









2 



.3. 5...(2k-1) 3. 5...(2k+1) 3. 5...(2k + 3)J 



and the equation (79) can be written 



i 

2 






2 



a? 



,6 



3.5 



I . * * 



X 



\ 



3~ 375 



l-f+--.., +2 1 






3 3.5 



3 



a? 






x 



o . D o . O . # 



/ 
/ -J 



= o, 



where x is written for sa. Now we have 



e**dx = e**lx-^+-^- 

\ t> O . 



and therefore 



9 4. 

«/y tAy 

¥ 377 



fa: 



— 1 — JU-2 
/v» »•/? a x 



<3 2 



Jkc* 



c/x 



Jo 
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and the equation may be written 



1-W1--J- 

x' 



x 2 



, 1 2 1 



1-2 

xe * 



e hx% dx = 





0. . 



(80) 



22. The left-hand member of equation (80), being equal to 

2 (l-fv) + 2(W) «*+..., 
is positive when # = 0. To determine its sign for large values of x we observe that 



•x 

e hx2 dx = 



rx 



e ix2 dx + x l .e ix<2 xdx 

Ji 



n 



rx 



e) x<z dx + x V^ 2 — #+ x 2 e w 'dx 

JO J ] 



fl 



•a; 



e ix2 dx + x V* 2 — e J +x d e hx2 — e* + 3 x V* 2 ote 

Ji 



. . . 5 



a; 



and therefore there is an asymptotic expansion^ for a?e * x2 e^dx when x is large in 
the form 

X "T" CC "T" OX ~T" O • 0£C "T" O . . / X "T"..». 

Hence the expression in the left-hand member of (80) is asymptotically equal to 



^(i-JA-L-i.-.Y,-*- 

V x ar \ r 



1 + 4+^7 + --^ + 






2T 



The term of highest degree independent of v is — 2^~ 4 ; the term of highest degree 
containing v is 8^~ 6 . It follows that the expression is always negative when x is 
sufficiently great. The expression therefore changes sign for some positive value of x, 
and the equation (80) has at least one positive root. 

23. When v = the equation (80) becomes 



1- 



x 2 -l 

X 



— l r 2 

e 2 



"X 

e ix2 dx = 0. 



If x 2 < 1 the left-hand member is necessarily positive. We shall take x 2 > 1 and 
write the equation 



x 



'X 



x 2 -l 



e^ 2 _ e* x2 dx = 0. 



Let y denote the left-hand member of this equation. Then we have 






O.v.2 



2x 






2 



e 



h* 



dx [x 2 —l x 2 —l (x 2 — l) 2 J (x 2 — 1)' 

Since this expression cannot vanish, the equation cannot have more than one positive 
root. 

* For the suggestion that this step might prove useful in demonstrating the existence of a root oi 
equation (80) I am indebted to Mr. G. H, HARDY, Fellow of Trinity College, Cambridge, 

2 c 2 
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Again, when v = \> the equation can be written 

"X 

x 2 e^ x2 dx \ = 0, 

*As \ tV tAj / |_ JO 

or 






O —1^-2 



\\ 6x l x e^ 2 + a?#*dx = 0, 
3« 4 -2a: 2 +l Jo 

where the left-hand member is certainly positive when x 2 < 1 ; also the differential 
coefficient of the left-hand member with respect to x is 

oX { oX — Zi ) 1^2 

and this expression cannot vanish for any value of x which is greater than unity. 
Hence the equation (80) cannot have more than one positive root. 
24. Now take v — 0, and write the equation (80) 

i / x 2 x 4 " \ 

— ( 1 — — + — — - — ...1 = 0. . . . . . . . (81) 



X JL \ o o * O 

When x 2 = 4, we have 

, x 2 , x A / 4 , 4 2 

1 — — + — — ... = J. — ■— - + 



o o . \ o o , 



44123 i' f ] _±\_ _ 4' / 4 



and 



o . . • » lo O . * . * X D \ J-// O . . , . 1 »7 \ £ J. 



1 1 45045 



X .1. o o , . . . X o 



Hence, when x 2 = 4, the sign of the left-hand member of the equation (81) is plus. 
When x 2 — 5, we have 

-r z — • • • =: 1 — "^ "T 



O o • o o » 



20056 15625 / 



3.7.9.11.13.3 3.7.9.11.13.3.17V 19 



Also 



3.7.9.11.13.3.17.19.21\ 23/ ""' 



15625 14 ^ „-„ 
> 677. 



17 19 



and therefore 



20733 



3.5 ° 3. 7. ..13. 3' 



but 



1(3. 7. ..13. 3) = 20270 + i < 20733. 

Hence, when x 2 = 5, the sign of the left-hand member of equation (81) is minus. 
It follows that the value of x 2 , or s 2 a\ which satisfies the equation is between 4 and 5, 
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25. Again, when v — \, equation (80) can be written 



3« 2 +l 


/\ 






' 1 - 


-l~ 


3a;*-2» 2 +l 


I 


3 3.5 



= 0. 



(82) 



If we put x 2 = 4, the left-hand member becomes 



41 V 3 3.5 



* • f • 



Now 



, 4 4 2 
1-- + 



645461 



8 



3 3.5 



+ 



and 



O • «-)•»• 1.0 O. O • • • J- / 

^(3. 5. ..15) = 642715+ — 
41 v 7 41 



4\ 

^lW" 1 "'''' 



Hence I — tt + o — -—...>—, and the sign of the left-hand member of (82) is 

O O • 4: 1 

minus when # 2 = 4. When we put x 2 = 3, the left-hand member of (82) becomes 

5 -d-U 32 



ii 



3 3.5 



but 



, 3 , 3 2 

1 — — + 



3 3.5 



3 9,3 
+ 



9 



and 



5 5.7 5.7 5.7. 11 

s L_fi_A+ ) 

7 5.7.11V 13 r 



11 7 V 3 3.5 



l ~h- 



or the sign of the left-hand member of (82) is plus when x 2 = 3 W It follows that the 
root of the equation (82) for x 2 lies between 3 and 4. 



Instability in respect of Displacements specified by Harmonics of the First Degree. 

26. When n > 0, we have to calculate expressions for A m B m C n , D n from the 
formulae of 8 17. If n = 1, we have 



s 2 a 2 



K = A, = 1 1 ( /l + ^ l) + a^ -F^-^- 2 (2E 1 -P' 1 + 2a 2 Q' 1 ). 
Now if we put h 2 = and P = 0, we find 



2.5 2.4.5 •"*• ' 2. 4. ..2k. 5"" 



# t/j — xjL 



3 . sW 



1- 

5 2 



7jV ,_ y +1 (2_K+3)^a^ 

2.4.5 '"^ ' 2. 4. ..2k. 5 



O V 



a' 



+ 



sV 



2.5 2.4.5 



-...(-)« 



s 2 K -2 a 2 K 



/j . 4 . . . iu/C . 



(83) 
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and therefore, if we put h 2 = and h 2 = in A 1} we get 



Xi-J — - 



\A/ ZlXX 


~ s 2 a? 
2 

*** 


2>v 5 




+ 2 A 
3i> 5 



s 4 a 4 
2T1 



— ££ _|_ 



•■• \ ) 



s 2x a 2 ' 



2. 4. ..2«: 



~2~ 



--(-)\: 



q2/c-2^2k 



= 0. 



2 . 4...(2k-2)*", 



It follows that Ax vanishes to the first order In h 2 and P, and therefore, as has been 
explained, we must evaluate the limit of AJi^ 2 when h 2 and F vanish. We have to 
expand the terms of f l9 s 2 1 and acW\/da correctly as far as h 2 ; in calculating the 
remaining terms of A 1? we may put h 2 and ¥ equal to zero in E l5 P\ and Q\. 

The terms of /J which are of the first order in h 2 are 



A/l 



a 



1 



7.9sV A 1\ , y 7.9. ..( 2k+5)s 2 «~~W" /I 1 , 

2.5 2.4.5.7\7 9/ *" 1 } 2 . 4. ..2k . 5 . 7..'.(2k+3)\7 9 2k+5, 



The terms of s 2 8 1 which are of the first order in h 2 are 



Alv 



2 



"3 1 



s 2 a? 



_0 DS Zi . 4 , 



, , v+1 7.9...(2k + 3)s 2 " 



2sy2lC 



2. 4. ..2k. 5.7.. (2k+1)\7 ' 9 



1 1 
+ r + ...+ 



2/c+3. 



Hence the terms of (/i + s 2 ^) a?/3v which are of first order in h 2 are 



It Lv 

>■— — ' I ■ ■ 

15i> 



,2^,4 



5* 



<r 9sV 

¥ + 27477 



... ( ; 



~2k— 2^ 2k 



< ju . 4 . o . Zi 



u 



K 



11 1 

7 9 2/c + 3, 



Again, when we keep those terms only which are of the first order in h 2 and F, we 
find from (52) of § 1 1, ' 

1 h 2 



p 2 = 



2 . 5s J 5^ 



2 e„2 ' 



F4=; 



^ 



2 . 4 . 5 5s 



2 ' 



P6 



^8 = 



Pl0 = 



2.4.6.5.7 



(A 2 +P), 



^J.4t.U.O»«J»# .«7 



7.9s 2 h 2 (± + h+F7s 2 



n . „ a * = = - - ( 7 . 9 . lls 4 /i 2 (^ + - + — ) +F7. 9.sH 
2.4.6.8. 10.5.7.9. 11 L \7 9 11 J I 



'± + ± + ±-\ 



* ? 



Pa< = (-)" ; 



2. 4. ..2«r. 5. 7...(2k+1) 



7.9...(2/c+l)s 2 *- 6 /i 2 :(i+I+...+ 1 



7 9 



2«r+l 



+ & 2 .7.9...(2k— l)s 2 



> ; 
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and therefore the terms of adP'Jda which are of the first order in h 2 and k 2 are 



_2_ 
3v 



s 2 A[2p 2 a 2 +4:p 4! a i +... + 2Kp 2K a 2K +...]. 



where p 2 ,... have the above values, that is to say, these terms are 



2Ah' 
15v 



a 



sW 



5s 2 2.4.7 



+ ...(-)' 



s^-V 



2.4...(2k-2)\7 9 



A l 

ZT + - + ...+ 



2/c+l 



+ 



2AF 
15y 



a' 



s 2 a° 



2.5 2 . 4 . 7 



+ ...(-)- 



^"V" 



"1 



2.4...(2k-2)(2#c+1)' 



It follows that the terms of the first order in h 2 and k 2 in 



are 



Aft 2 a 



2^2 



3i/ 



L_ ^L + 



A 4 



15i/ U* ■ 2 2.4 



<ia 



—•(-)* 



s 2 *~ 2 a 2 * 



2. 4. ...2k 



or they are 



2„2 



+ 



2AF« 
15v 



L5 2 7-2. 4 + "^ ; 



s 2k - 2 « 2k 



2k+3 2.4...2k 



APa 2 _ 
15s 2 



£sW 



+ 



2 AFa 2 /l s 2 a 2 1 s 4 a 



4™4 



15A \5 2 7 2.4 



+ ... 



Again, when we put h 2 and k 2 equal to zero, we find 



2E 1 = 



-Fi = 



2a 2 Q' 1 = 



2A 
5s 2 

2As 2 
15i/ 

2 As 2 



15v 



1 + 



2? 



s 2 a 2 



s 4 a 4 
271 



+ ••• (~) K+1 



s 2 "a 2,c 



2. 4. ..2k 



a' 



+ 



s 2 a 6 



2 .4 2.4.6 



-•••(-r 1 



s 2 '- 4 a 2 " 
2. 4. ..2a:' 



a 2 -la 
+ 



2„6 



1 s z a' 



5s 2 7 2 9 2.4 



...(-)« 



6 >-4 a 2 K 



(2«t+3)2.4...(2k-2)"" 



and therefore the terms of the first order in h 2 and k 2 in 



2„2„2 



k 2 s 2 a 
3s 2 +h 2 



(2E 1 -P' 1 + 2a 2 Q / 1 ) 



are 



1 
3 



k 2 a 2 



2A 
5s s 



e 



k'o> 



+ 



2A f s 2 a 



sV 



g* K a 2K 



5sVlO~~2.4.7 '"^ ' 2T4lT72^c+3) 



Hence the terms of A x which are of the first order in A 2 and A 2 are 

1 f^fc ^ e ~\&(& 



and 



1 • x\_i J\.Cf/ 

iim A = T7 = — r- e 
nr 5 s rp 



•&**«* 
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27. Again, when we put h 2 = and ¥ = 0, we find that 



Bx = 



Sp a da da 



where 



i(Ms%) 



2A 
I5u 



s 2 a 2 



!_— + 



s*a* 

271 



-(-) K : 



s^a 2 " 



2. 4. ..2k 



2 dgi 2A 



s 2 a 2 



+ 



s 4 a* 



2 2.4 



— •(")" 



^2/c a 2« 



2. 4. ..2k 



a ^Qi-Q' 1 = A 



cZa 



15^ 



1 , sW 
+ 



,2k 



5 '2.7 2.4.9 '" V ' 2.4...2k.(2k+5) 



_A_ 

15v 



'6_iU^-7£-iW^YS-i 



2 \7 



2 V9 



(")* 



s 2K a 2K 



6 



2. 4. ..2k \2k+5 



1 . 



A _>*,„ . 2A /l 1 s 2 a 2 . 1 sV 



15v 



e-**" + — ( 



+ 



5v \5 7 2 9 2.4 



Hence we have, when h 2 = and A 2 = 0, 



1 = - — ( 1 — — ) e 2 + 



< 



2v, 



5v 



l_lsV l sV 
5 7 2 92.4 



-t 12 



2A 



•sa 



J> 



5p K } Jo 



(85) 



Also when w = 1 we have 



and therefore 



and we have also 



C w — G r — g ^n(^): 



lim 



Ox 



a' 



* =0 /^ 2 



3s 



2' 



D„ = D x = F [^(Aa) +3^ a (*a) + (1 h-A 8 /* 8 )- 1 ^^)], 



(86) 



and therefore 



r D x 3/-l\ , 1 



3. 5 



10 



(87) 
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28. In the case where n = 1 the condition of gravitational instability, viz., 
lim k=0 (A 1 h'' 2 I) 1 k'' 2 — Bidh^k' 2 ) = 0, becomes 



O J\.d __!£.€ 



10 bs 2 v 



e 



£S%2 



+ 



a 



^ij( 2 -~y^^ 



sa 



0. 



(88) 



But we have 



*sa 







x'e'^dx = -(A 3 +3m)e' |A2 + ZeT^dx, 

Jo 



'sa 



and therefore the condition of gravitational instability becomes 



sa 

e-Wdx-e-***} Ssa+ (saf + 



19 



l 



L • — • \20 V )W\ 



0. 



If now we put 
the equation becomes 



s 2 a 2 






2 LV 



2 



e~* 2 c& 



77 Jo 



— e 



- Z 2 



19 



32 + 22 3 +(-^--4*/)2 5 



= 0. 



. (89) 



where the factor 2tt l has been inserted because the expression in the square brackets 
is tabulated in many easily accessible books. 

Let y denote the left-hand member of the equation (89). When % is small, y is 
small of order z b . In fact, we have 



e~ p dt = e~ z " 







. 22 3 , iz b s 
o o . 5 



k^2k+1 



2*2 



3. 5...(2/c+l) 



+ ... 



and when 2 is small, the first approximation to y is -(3-4i>)z 6 . When z is great, 

e~ p dt is approximately equal to |vV; and thus y is positive when z is great. The 

equation (89) has a root zero and at least one positive root. The zero root is 
irrelevant to our problem; it is introduced in transforming equation (88) into 
equation (89). Now we have 



~f = -*V* 2 j(l5-20z,)-|(l9-20z/)2 



2 



and, since 15 — 20^ and 19 — 20i> are positive when y<f, the expression last written 
vanishes for one positive value of z. Hence it follows that the equation (89) has only 
one positive root, and there is one and only one positive value of s 2 a 2 which satisfies 
equation (88). 

By means of the tables it can be shown that, when v = 0, the root z lies between 
1'9 and 2, so that s 2 a 2 lies between 7*22 and 8. When v = !, the root z lies between 
1-8 and 1-9, so that s 2 a 2 lies between 6*48 and 7*22. 

vol. ccvu. — a. 2 D 
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Stability in Respect of Displacements Specified by Harmonics of the Second and 

Third Degrees. 

29. When n — 2 and h 2 = 0, W — 0, we have 

A n = A 2 = f!(/ 2+s 2 2 ) + a ^ 2 + f 2 +3E 2 +3a 2 Q' 2 , 

5v da 5 

where 



J 2 — A 



1 



8 



8.10 



<A, 2 4- "•- LW *W _ / _ V 8.10...(2K+6)s ac a 2 « 
2.7 2.4.7.9 "^ ' 2.4...2ic.7.9...(2k+5)'". 



0, = A 



5,i2 8s 2 a 4 / v+1 8 . 10...(2«:+4)s 



2k-2q2k 



6s 



2.4.7 



2.4...2k.7.9...(2#c+3)' 



9. = A 



'2 



1 of 



,2 K -2 a 2K 



___ 8s 2 a 4 (_y+i 8 .10...(2k+4) s 

6s 2 ' 2.7 2.4.7.9 "^ ' 2T4...2*. 7.9 ...(2k+5) **'J' 



p/ 



Q', 



2A^ 
5t> 



a' 



+ 



a" 



4 2K~4 r .2/c 



8sV , y 8.10...(2ic + 2)s 

12s 2 ' 2.4.7 2.4.6.7.9 *"'* ' 2.4...2ic.7.9...(2ic+3)*"_ 






a' 



>2K-2~2ie 



8s V / y +1 8.10...(2k + 4)s 

6. 7. s 2 ' 2.7.9" 2747779T11 "^ ^ 2.4...2k.77977.(2k+7)"V 



From these we find 



Ao — — 



-2 



2Aa 
5v 



2 r 



1 aV &sV , y (2/c+2)(2/c+4)s 2K a 2 K 

L6 2.7 2.4.7.9 '"^ ' 2.4.6.7. 9... (2k+5) 



1 



sV 



+ 



6.7 2.7.9 2 



8s 4 a 4 / y (2K + 2)(2K+4)g 2 "a 2K 

.4.7.9.11 '"^ ' 2.4.6.7.9...(2k + 7) 



3A 



2A 



"1 

6 

w 

"1 
6 



~2 2 



+ 



8s*a 



4^4 



2.7 2.4.7.9 



/ y (2* + 2)(2;c + 4)s 2 «a 2 « 
"^ ' 2.4.6.7.9...(2ic+5)"'. 



s 2 a 2 



+ 



8s*d 



4^4 



2.5 2.4.5.7 



■MM* f IWW W W W j 

••• V / 



, (2«r+2)(2<c+4)g 3K a 2,: 
2.4.6. 5.7... (2k+3)'"_' 



By means of the identities 



(2k+2)(2k+4) 


(2/c + 3)(2k+5) 
(2k+2)(2k + 4) 


(2k+5)(2k + 7) 
(2k+2)(2k+4) 


(2k + 1)(2k+3) 



1- 



+ 



2k+3 (2/c+3)(2k+5)' 



= i 



6 






15 



2«r+5 (2k + 5)(2/c + 7)' 



= 1 + 



1 



2k+1 (2k+1)(2k+3)' 
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we transform the series 



1 
5 



1 
6 



x' 



+ 



8x 4 



1 

5 



2.7 2.4.7.9 

.2 



( - V (2« + 2)(2>c + 4) 



x 



2k 



2.4.6.7.9...(2/c + 5)'"J' 



x 



+ 



6.7 2.7.9 2 



8 a; 4 , y (2k + 2){2k + 4:)x 2k 

.4.7.9.11 '"^ ' 2.4.6.7.9...(2/c + 7 



(2k + 7) J 



1 
6 



a? 



+ 



8*' 



2.5 2.4.5.7 



-...(- V J2k+2)_(2k + 4) 



a; 



2k 



2.4.6.5.7...(2/c + 3) J* 



respectively, into the forms 



1 



16 
16 



2 

3 



1_± 



1 



3.5. 



•a^fl- 2 + 



6 






15 



3.5 3.5.7, 



— x 



3.5 3.5.7, 
J 1 6 



V3.5 



+ x" 



+ 



1 



+ 



3.5 3.5.7 3.5.7.9, 
15 



3.5.7 3.5.7.9 
1 6 



+ x 



.4 



+ 



15 



3.5.7 3.5.7.9 3.5.7.9.11 



16 



1 + 2 



1 



x 2 1 + 



3 3.5 



• 4 1 I , 2 



1 



3 3.5 3.5.7 



or 

J_ 
16 

J_ 
16 



2 4 



/I o^ 2 T 4 

2(- — — - 1- 

3 3.5 3.5.7 



-... +3 



x 



+ 



CC 



203 



\3.5 3.5.7 3.5.7.9 



1 



x 



2 



+ 



X 



3.5 3.5.7 



6 



x 2 



+ 



# 



4 



3.5 3.5.7 3.5.7.9 

+ 15' ] 



a? 



+ 



a; 11 



3.5.7 3.5.7.9 3.5.7.9.11 



1 



16 



1- 



4 6 

ar + — — 



^ 



XT 



,x 4 



rt +...+21 + 

3 3.5 / \ 3 3.5 



a 



X 4 



+ 



X 



,3 3.5 3.5.7 



Now we have 



iAj %Aj 

~3 + 3T5 



1 
3 



ar 



+ 



x- 



3.5 3.5.7 



•a; 



x x e ^ \ e ix2 dx, 





*x 



x* 



•x 8 e lx<2 1 e hx<i dx, 





and therefore the three series are respectively equal to 



16 



1 + —2 + -7 ) X le 

iAj vis 



■w 



X 



e^dx— ( -9+ -i 
o \ar ar 



16 



1 , 15 

—~A + 



■a; 



X 



X 



6 



k ar x 4 ar/ Jo 



16 



*x 



1- 



cc 2 — 2 ^ja; *e hx \ e* x2 dx 



x" 



x 
2 d 2 



o 
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On substituting for the three series in the expression for A 3 we find 



A< 



Aa 2 


" 1 


8v 


_s 2 a 2 


A 




16s 2 



45 

s Q a e 



1- 



15 45 \ 



s 2 a 2 



sW 



s 6 a Q 



. ) (sa) 1 e h 



s 2 a2 



s 2 a 2 «V 



32 45 



sV 



sV 



"sa 

e^dx 











(m)-^-**" 


m sa 

e ix<2 dx 




%/ 






• (90) 



30. Again we have, when h 2 = and F = 0, 



B 



2 






8 <ZE, + dtf,_Q, 



5ct c£a 



da 



2 ) 



and with the values already used for f 2 , . . . this gives 



R 



A 
5v 



s 2 a 2 



+ 



8s\i 



4 ,4 



6 2.7 2.4.7.9 



-...(-) 



, (2k+ 2) (2 K +i)s 2K a 2K 
2.4.~6T7. 9...(2/c+5)""_ 



+ 



8A 
5 



1 



8s 2 a 



2„2 



4„,4 



+ 



8. 10s 4 a 



7 2.7.9 2 



.4.7.9.11 v ; 



« (2/c+2)( 2 /c + 4) (2k+6)s 2ii ct 
2.4~.6. 7. 9~...(2k + 7) 



~5> 



s 2 a 2 



3.8sV _ , _ y (2fc-l)(2fc+2)(2K+ 4)g 2 'a 
6.7 2.7.9 ' 2.4.7.9.11 '" ^ ' 2 . 4 . 6 . 7 . 9 ...(2k+7) 



2k 



+ 



The first of these series has already been transformed into 



16v 



1 + -s-5 + -T-: ) (*«) l e lsW I e^ 2 £&■ 



u \ 



s 2 a 2 



sV 



pa 

Jo 



1 



3 



2 2 ' 4 4 



By means of the identity 



(2k+2)(2k + 4)(2k+6) _ 1 3 



+ 



9 



15 



(2k+3)(2k+5)(2k+7) 2k + 3 (2k+3) (2k+5) (2k+3) (2/c+5) (2k+7) 



we transform the series 



into 



5 



8x 



2 



7 2.7.9 



+ -(-)* 



(2kt+2)(2k + 4: )x 2 ' 
2.4.6.7. 9..."(2k+7)"'_ 



16 



2 4 

3 + 3T5 



3 



x' 



+ 



X 



3.5 3.5.7 



+ 9 



1 



x 



+ 



X 



o.o o . o . / o.o./.y 



15 



ar 



+ 



03* 



3.5.7 3.5.7.9 3.5.7.9.11 



which is the same as 



16 



.,3,9, 15\ - x 



'a; 







1 



+ —T + 



15' 






# 



6 
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and by means of the identity 

(2k-1)(,2*+2)(2k+4) _ x 
(2k + 3)(2k+5J(2k+7) 

we transform the series 



10 



w* - * 



51 



120 



2k + 3 (2k+3)(2k+5) (2k+3) (2k+5) (2k+7) 



1 
5 



X' 



3 . 8x* , u( 2K-l)(2^+2)(2/c+4)a; 2 ' t ' 

6.7 2.7.9 2.4.7.9.11 '"^ } 2 . 4 . 6 . 7 . 9 ...(2k+7) '". 



into 



16 



+ 



o . D 



— ... 1—10 



iA> «A> 



O O , D O • O . I 



-...1 + 51 



1 



/y>^ 



+ 



af 



o . D o . . # o . • / • «/ 



120 






: + 






3.5.7 3.5.7.9 3.5.7.9.11 



which is the same as 



16 



10 . 51 . 120 



1+^+— 7 + 






A 



X 



6 



X 



■i 6 -** jV cfa- f 1 + 1} + 12 ° X 






Hence we have 



R,= 



16i> 

A 

2 



3 



25 . 240 



s 2 a 2 



"I r~l + 



s 6 a 6 



3 + 



22 . 105 . 240\ 



2 2 

O W 



J. 



s 4 a 4 






s fl a fl 



)(sa) l 



e 



■IsW 



sa 

e^ xi dx 



o 



i + JL + 44_(' 1+ JL, + J, + ^( M )-.^™ 



s 2 a 2 



s 4 a 4 



$W 



sV 



s 4 a 4 



s 6 a 6 



e ia:9 da? 



• • 



(91) 



Affain we find 



li'm Zl — x 

J=, F" 5.7.9.S 2 



(8s 2 a 2 +175t>), . 



» • « • » 



anc 



lim A 



D 



2 



33 



P ^ 7 Q 



• •••*•• 



• • 



(92) 



(93) 



Hence the equation lim A=0 (A 2 D 2 & 3 — B 2 C 2 ^ 2 ) = is 



(8x 2 +175v) 



3 , 25 . 240 1 



+ -r + 






— 8v 



4 



8+22 + 105 + 240)*-^ 



x- 



x 2 ' a; 4 * # 6 / ^o 






cc 2 a? 4 a; 6 / \ x 2 x* x % 



'X 

6^X fj/yt 



Jo 



+ 33 



2 



9 A)„( 2 x 2 -2^ 
xV \ x 2 



90' 

x 4 / 



% x 



-l*-h& 



M* 



x x e **" ( e 2 ""~ <ia? 





V i 



2 



17_45' 

eA> fA> 



o 2 1 A Y ~"1 



•^ 



e^dx 



o 



= 0, . (94) 



where # is written for sa. A factor *ry- 3 has been omitted in forming the equation, 
as neither v nor s is supposed to vanish, The terms of lowest degree in the left-hand 
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member of the equation (94) will be found to be ^80^(24^—19), which is negative 
when § > v > 0. Hence the left-hand member of this equation is negative when 
x — 0. Also it will be found by the method of asymptotic expansion (of. § 22 above) 
that, when x is great, the left-hand member is approximately equal to — 1280#r 4 for 
all values of v. After the previous cases, in which the corresponding equations have 
a single root, we are led to expect that in this case there is no root, for it is unlikely 
that there is more than one. We proceed to verify this expectation in the cases 



where v = and v =- l 



4* 



31. Multiply the left-hand member of (94) by x* and put v = 0. We get 



8 



(3^ 4 +25^ 2 +240)---(3x 6 + 22^ 4 +105x 2 +240)ir" 1 e"^ 2 e**dx 



X 



+ 33 



"x 



(2« 4 -90) - (2« 6 -2x 4 -30,x 2 -90) x~ x e~^ e ix2 dx 



0, 



or, since 10 is a factor. 



•x 



9* 4 +20x 2 -105-(9a; 6 +ll£c 4 -15x 2 -105)«~V i:c2 & x "dx = 0. 

Jo 

The term of lowest degree in the left-hand member, when expanded in powers of 
x, is — - 3 % 6 ; when x is great, the left-hand member approximates to —128. Now 
multiply by xe ix \ and put 

rx 

y = (9£ 5 +20x 3 -105«)e^-(9* 6 +ll« 4 -15«: 3 -105) e&dx. 



We know that when x is small y is small and negative of the order — ^ic 7 , and that 
when x is great y is great and negative of the order — 128,*e Jj;2 . Now 

^L = (5ix i -S0x 2 )e iz '-(5Ax 6 +Ux 3 -S0x) V ePdx, 
dx Jo 

and, if we put % for x~ l dyjdx, 

z = (54x 3 -30^)e^ 2 -(54^ 4 +44x 2 -30) ['e&dx, 

Jo 

where z is negative both when x is small and when x is great ; also 



dz 
dx 



rx 



88^V* 2 -(216x 3 +88;£) e*"dx, 

Jo 



and if we put w for x x dzjdx, 

w = 88a^ 2 - (216^+88) X & x "dx, 

Jo 

where w is negative both when x is small and when x is great ; and now 



dw 

dx 



= -128are**-432a5 



e* x2 dx. 



Jo 
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which is always negative. Hence w is always negative, and therefore dzjdx is always 
negative, therefore also z is always negative and dyjdx is always negative, and 
therefore y is always negative. Thus the equation y = has no real root other than 
the irrelevant root x = 0, which was introduced in the process of forming the 
equation. 

Again, when v = \, we multiply the left-hand member of (94) by x Q and obtain the 
equation 

(23x 6 +128^-70^ 2 +3675)-(23cc 8 +105^ [* e*"dx = 0, 

Jo 

of which the left-hand member is of the order —x Q when x is small and —x 2 when x 
is great. We put 

y = (23x 7 +128£ 5 -70x 3 +3675^)#^ 

Jo 

and then we put 

__ 1 dy _ 1 dz __ 1 dw 

%Kj \AitAj *aj \Aj%Aj t/y \AitAj 

and find 

~ = -(512cc 4 +4552cc 2 + 3640)e^-8832 ['<*"dx, 

ax j o 

which is always negative. Just as before, we deduce that y is always negative. 
It is therefore proved that the equation 

lim, =0 (A 2 D 2 r 2 -B 2 C 2 ^" 2 ) = 
has no real root. 

32. When n = 3, and h 2 = and k 2 = 0, we have 



A 3 = ~(Ms%) + fa + a **JL + |(4E 3 +4a 2 Q' 3 ), 



where 



s 2 a 2 . sV , x „ s 2 *a? K 



/, = A l-^L + ^.-...(_)« 



2.4 v ' 2.4...2«r 



Q -A/, s 2 a 2 , sW , u s 2 *a 2 * 



s 2 V 2 ' 2.4 """ v ' 2.4...2/c'*'" 
F= 4A/ a 2 . sV sV . / ^ s 2x ~ 2 a 2x 



Iv 



7" 279~2.4.11 ^"'' 2.4...2k.(2k+5)"' 



i^3= ^1 -= + 



7 2.9 2.4.11 + " -1 j 2.4...2k(2k + 7)"7' 



Q' 3 = 3 M- J- + ^ **_ + ...(-)-« 



7v\ 7.9 2.9.11 2.4.11.13 '*' v 7 2.4...2k(2k+5)(2k+7) 
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We use the identity 



2 



(2k + 5)(2k+7) 2k+5 2k+7 



to transform Q' 3 and then we use the results 



We thus find 



1 — ■ h 

2 2.4 



e -isW 



s 2 a 2 



+ 



sW 



*sa 



7 2.9 2 4. 11 



— ... = (sa) 7 \ x 6 e * x2 dx, 



s 2 a 2 



+ 



4 4 



9 2.11 2.4.13 



— . . . = (sa) 



-9 







sa 






A»= - 



-3 



4A 
7^ 2 



*s« 



3(sa) 5 A ^ 2 cfcc — 2(sa) 7 as 8 e ^(ix 

Jo Jo 



sa 



-{•ve 



. 28 / v 



*5Ct 






• (95) 



Again we have, when h 2 = and F = 0, 



B 3 = -i-(/ 3 +^3) + ^-^ 3 +«^- 3 -Q' 3 , 

7v 7 a da da 



and here we have 



a 



d 3!i 

da 



Q' 8 =- 



3A 



7v L 7.9 2 



1 ^ <x ds a / v i-^k — i 

.9.11 2.4.11.13~'"^~' 2.4.. .2^2^ 



(2«:-l)s 2 'a 2K 



(2k+7)(2k + 9)'*'. 



while the other series can be obtained from those written above, 
identity 



We use the 



2«r-l 



(2k+7)(2/c+9) 



+ 



2k+7 2k+9 



to transform the series last written, and we use also the results which we used in 
obtaining the expression (95) for A 3 , and find 



B a = 



A 
7v|_ 



12 (sa)" 



rsa 



6,» l&dx. 



x"e 



o 



sa 



15 (sa)~ 9 x 8 e-* x2 dx+l0v(sa) 

Jo 



9 



sa 



x 8 e » x<i dx 







• (96) 



We find also the results 



1" ^3 



o,o, § . y . s 



(2sV+49i>), 



* • e • » 



iim* ~. a , - — — — 

k 2 4.5.7 



• (97) 



(98) 
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Hence the equation 



Km A=0 (A 3 D 3 &- 2 -B 3 Cyr 2 ) = 



becomes 



3 



3x 2 -\ v\x 7 x 6 e hx2 dx-\-ve * : 



o 







2x 7 x s e lx<2 dx 



- -L (2x 2 + 4:9v) 

27 v 



12.x' 



a; 



x 6 e ix<2 dx—(l5 — l0v)x 9 cc 8 e ^cfe 

. .0 'Jo 



x 



0, . (99) 



where a? is written for sa. 

33. An irrelevant factor v has been introduced into the left-hand member of (99), 
We find that when x = this expression becomes 

^539-^ 



81\ 



3 



which is positive for all admissible values of v. We find also, by the m#fchod of 
asymptotic expansion, that the expression is positive, when x is great, for all values 
of v. We proceed to show that, in the important cases v = and v = J, the equation 
has no real root. The left-hand member of (99) being an even function of x> we may 
treat x as positive* 

When v as 0, the left-hand member of equation (99) is 



9 



x 



f' 

Jo 



7x 2 aJV** dx-2 



X 






H 



which is positive for small values of x. The differential coefficient of the expression 
within the square brackets is 



I Ax 



•x 

x 6 e~' ix2 dx + 5x s e~^ x \ 



Jo 



which is positive for all real values of x. Hence the left-hand member of 
equation (99) cannot be negative if x is positive, or the equation has no real root. 
When v = J, the left-hand member of equation (99) becomes 



x 



•9 



'49 
V9 



x 



133 a 

- — -X 

9 



rx 



x 6 e sx2 dx 



/.0 



'62 2 

— X — 

27 



1225 \ 
54 /. 



x s e ix2 dx 







+ - er** 
4 



0. (100) 



The expression within the square brackets is greater than 



5 (« 4 -3cc 2 ) \ X x 6 e- ix2 dx-(Sx 2 -22) \"a^e~^dx, 

Jo J o 



. . (101) 



where this expression is obtained from the other by replacing every positive coefficient 
by the next smaller integer and every negative coefficient by the next greater 
integer. 
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Since 

'x rx 

x s e~~* x2 dx = —x 1 e" ix ' 2 +7 x % e~~ hx<2 dx, 

o Jo 

the expression (101), when multiplied by 7, is 

rx 

(5^ 4 -36^ 2 +154) x*e-** dx+5(a?-S)x*e-** (102) 

Jo 

Also we have 

* /I /y.2 ^4 \ 

o^ ^ = " e U + ^TI + 9711713 + •••]'* 

and therefore the coefficient of x 9 e~ ix2 in the expansion of (102) is -^(154 — 9 x 15), the 
coefficient of x ll e~^ is ~ 9 - 9 -(154 — 36 x 11 + 5 x 99), and the coefficient of x^ +2K e~ hx2 for 
all values of k which are greater than 1 is 

154 36 , r 

+ 5 



m 9.11...(2*+5)\(2fc + 7)(2#e+9) 2k + 7 

or 

20k 2 +88/c+145 

9. ll...(2ic+5)' 

Hence all the coefficients are positive, the expression (102) is positive, and the left- 
hand member of (100) is positive for all positive values of x. Thus, in this case also, 
the equation (99) has no real root. 

Summary of the Solution of the Mathematical Problem. 

34. We have now solved in essentials the mathematical problem of the vibra- 
tions of a gravitating sphere, initially homogeneous and in a state of hydrostatic 
pressure, and have found the conditions of gravitational instability. We have shown 
that, when any normal, or principal, vibration is taking place, the dilatation at 
a distance r from the centre is specified by the product of a certain function of r and 
a spherical harmonic of positive integral degree. We have shown further that, in 
each . such mode of vibration, the components of displacement can be expressed in 
terms of the same spherical harmonic, and that the radial displacement at a point 
distant r from the centre is the product of a function of r and the same harmonic. 
We obtained the form of the frequency equation, and the forms of all the functions 
which enter into its expression. 

We proceeded to investigate the conditions which must hold in order that the 
frequency equation may be satisfied by a zero value of the frequency. We showed 
that, when such a value is not introduced irrelevantly in the process of forming the 
equation, its occurrence points to genuine gravitational instability. We found that 
the condition of such instability is the condition that a certain equation, containing 
the variable quantity s 2 a 2 , or |7ryp 2 /(^ + 2/x), may be satisfied by a real positive value 
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of this quantity. The constant fi denotes the rigidity and X + f/x the modulus of 
compression. When the harmonic specifying the vibrations is of zero degree, that 
is to say, when the vibrations are radial, we found that the critical value of s 2 a 2 
lies between 4 and 5 if v, or /x/(X+2/x), is zero, and that it lies between 3 and 4 if v 
is \. In the case of vibrations specified by spherical harmonics of the first degree, we 
found that the critical value of s 2 a 2 lies between 7 '"22 and 8 if v = 0, and it lies 
between 6*48 and 7*22, if v = \. In the cases of vibrations specified by spherical 
harmonics of the second and third degrees we found that there is no critical value of 
s 2 a 2 , or that the sphere is stable, in respect of the corresponding types of displacement, 
for all values of X + 2/x. It was to be expected that the critical values of s 2 a 2 would 
increase rapidly as the complexity of the type of vibration, specified by the degree of 
the appropriate harmonic, increases ; and we appear to be justified in concluding that 
instability cannot occur in respect of displacements specified by spherical harmonics of 
any degree higher than the first. 

35. The result that the critical values of s 2 a 2 are lower when v = | than when v = 
means that a higher value of the constant X-f-2/x is required, to secure stability, 
when there is considerable rigidity than when there is very little rigidity. This 
result accords with general dynamical principles ; for it is a known result, and one 
which has been shown to be in accordance with such principles, that the frequency of 
any mode of vibration, involving compression, of a sphere free from gravitation 
diminishes as (X-f-2/x)///, diminishes, that is, as v increases. # Consequently, for a given 
value of yp 2 a 2 > the value of yp 2 a 2 /(\ + 2/x) which would be required, in order to reduce 
the frequency to zero, diminishes as v increases, or the critical value of X+2/x increases 
as v increases. 

36. The result that the critical value of s 2 a 2 is lower when n = than when n = 1 
means that a higher value of X+2/x would be required, to secure stability, in the case 
of radial displacements than in the case of displacements specified by spherical 
harmonics of the first degree. The spherical body of uniform density could be stable 
in respect of all types of displacement except radial displacements. If the value of 
s 2 a 2 were intermediate between the critical values corresponding to n = and n = 1, 
this would be the case ; and the body would tend to take up a different configuration, 
in which the density would be more concentrated towards the centre. The result 
that, in the case where n = 1 also, there exists a critical value of s 2 a 2 , which 
is not more than twice as great as the value associated with n = 0, the initial 
state in both cases being one of uniform density, suggests very strongly that there 
would be a critical value of X + 2/x, in respect of the case n = 1, even if the configuration 
were such that the body was stable as regards radial displacements. We should then 
have a body with a spherically symmetrical distribution of density, but with elasticity 
too small for this configuration to be stable in respect of displacements specified by 
spherical harmonics of the first degree ; and it may be inferred that the critical mean 

* Of. H. Lamb, loc. cit, ante, p. 173. 
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value of X + 2/x for such a body would not be very different from the critical value 
obtained for A. + 2/x by treating the body as homogeneous, and paying attention to 
those types of displacement only which are specified by spherical harmonics of the 
first degree. 

37. If this conclusion is admitted, as I think it must be, it would follow that a 
spherical planet with a spherically symmetrical distribution of density, and stable as 
regards radial displacements, might be unstable as regards displacements of the type 
in question ; and then it would tend to be displaced in such a way that the boundary, 
or any concentric sphere, moves to a position in which its centre no longer coincides 
with the centre of gravity, while the matter in a thin spherical layer becomes 
condensed in one hemisphere and rarefied in the other. The density being in excess 
in one hemisphere and in defect in the other, and the excess or defect at any point, 
at a stated distance from the centre, being proportional to the distance of the point 
from the bounding plane of the two hemispheres, the distribution of density may be 
aptly described as " hemispherical," and the state of the body may be described as one 
of "lateral disturbance." The concentration of density towards one radius, on which 
the centre of gravity lies, has the effect of diminishing the potential energy of 
gravitation, and this diminution may more than counterbalance the increment of 
potential energy due to strain. The proved existence of a critical value for X + 2/x (in 
the case of a homogeneous body) indicates that this state of things really can occur. 
An illustration of the nature of a hemispherical distribution of density will be found 
in §§ 47, 48 below. 

38. The results found by Jeans (1903) in the solution of the problem of the 
gravitating sphere subjected to an external field of force, which balances gravitation 
throughout the sphere when it is at rest, may be compared with those obtained above 
in the case where the gravitation is balanced by initial pressure. In Jeans' solution, 
just as here, the modes of vibration are specified by the spherical harmonics which 
enter into the expression for the dilatation ; and, in any normal mode, the formula 
for the dilatation contains a single spherical harmonic, and the radial displacement 
at any stated distance from the centre is proportional to the same harmonic. If the 
degree of the harmonic exceeds zero, instability can occur for a sufficiently small 
value of the resistance to compression, whatever the degree of the harmonic may be. 
It is not restricted to the case where the degree is unity, as it is in our problem of 
initial stress ; but the value of the resistance to compression required for instability 
diminishes rapidly as the degree of the harmonic increases. Instability enters first 
when the harmonic is of the first degree, # that is to say, for lateral disturbances. 
The critical values of s 2 a 2 are 6'72 when v = and 5*33 when v = \, the degree of 
the harmonic being unity. Since these values are a little less than the critical values 
found in the solution of the problem of initial stress, it may be concluded that the 
effect of initial stress, as compared with that of an external field of force, is to 

* The question of radial instability, was not considered by Jeans, 
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increase slightly the stability of the body in respect of disturbances specified by 
harmonics of the first degree, and to increase it enormously in respect of disturbances 
specified by harmonics of higher degrees, 



Application to the Problem of the Gravitational Stability of the Earth. 

39. For a body of the same size and mass as the Earth, the values of a and p in 
C.G.S. units are 6*37 xlO 8 and 5*53; the value of y being 6*65 x 10~ 8 , the value of 
\iryp 2 d 2 is 3*46 xlO 12 . In the following table the first column gives a value of s 2 a 2 , 
the second column gives the corresponding value of X + 2/x (the body being of the same 
size and mass as the Earth), the third and fourth columns give the values of the 
corresponding moduluses of compression in the cases where v = and v = \ , irrelevant 
entries being omitted. These moduluses are denoted by k and k lt The quantities 
given in the fifth, sixth, and seventh columns are the moduluses of compression of 
steel, glass, and mercury (denoted by k s , h g9 Jc m ). 
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According to these results, a homogeneous solid body of the same size and mass as 
the Earth, with a modulus of compression as great as that of steel, would have 
complete gravitational stability. If the modulus of compression were equal to, or less 
than that of glass, the planet would be unstable as regards radial disturbances, and a 
concentration of density towards the centre would take place, If the critical value of 
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X + 2/x, which was found in the case of lateral disturbances, is assumed to be the 
critical mean value of X + 2//, for a planet in which the mass is condensed towards the 
centre, then we may say that, if the mean modulus of compression were about equal 
to that of glass, and there were very little rigidity, the planet would be unstable as 
regards lateral disturbances ; but, if there were considerable rigidity, it would be 
stable. If, on the other hand, the mean modulus of compression were decidedly less 
than that of glass, though not so small as that of mercury, the planet would be 
unstable as regards lateral disturbances, even though it possessed a considerable mean 
rigidity. 

40. In order to settle the question of the gravitational stability or instability 
of the Earth, we must assign the appropriate values to the constants X and /x. 
Lord Kelvin's theory of elastic tides in a solid sphere led to the result that the tidal 
effective rigidity of the Earth is not less than that of steel. This result suggests 
that jjl should not be taken to be less than 8 '19 x 10 11 C.G.S. units ; but, since it was 
obtained by treating the Earth as incompressible, it affords no means of determining 
the value of X. Jeans (1903) proposed to deduce the values of X and ju, from the 
observed velocities of propagation of earthquake shocks. In a homogeneous elastic 
solid body, free from gravitation and initial stress, irrotational waves of dilatation are 
propagated with the velocity [(X + 2/x)/p ] i , where p is the density, and equivoluminal 
waves of distortion are propagated with the velocity [ji/pof, while waves of a third 
type are propagated over the surface with a velocity approximately equal to 
(0*9) [fi/pof- When a great earthquake takes place, the disturbance received at a 
distance from the source consists of three sets of disturbances : two sets of 
" preliminary tremors," and the " main shock." The first set of preliminary tremors 
is received at distant places at such times as it would be if it travelled directly through 
the Earth with a velocity of about 10 kiloms. per second. The second set of tremors 
is propagated apparently in a rather less regular fashion, but the times at which it 
can be observed at distant stations are nearly the same as they would be if it travelled 
directly through the Earth with a velocity of about 5 kiloms. per second. The main 
shock is received at distant places at such times as it would be if it travelled over the 
surface of the Earth with a velocity of about 3 kiloms. per second. # The identification 
of the three sets of disturbances with the three sets of waves which are theoretically 
known seems to be inevitable, and the discrepancy between the ratio of velocities of 
equivoluminal and superficial waves and the ratio of velocities of the second set of 
tremors and the main shock may be explained by the supposition that, while the 
velocity of transmission of these tremors depends upon the mean rigidity of the Earth 
as a whole, the velocity of transmission of the main shock depends upon the average 

* Eeference may be made to a Memoir by R. D. Oldham, " On the Propagation of Earthquake Motion 
to great distances/' London, 'Phil. Trans. Roy, Soc./ ser. A, 194, 1900, and to the Reports of the 
Seismological Committee of the British Association, in particular that published, in 'Brit, Assoc, Rep,/ 
1902, 



GRAVITATIONAL STABILITY OF THE EARTH. 215 

rigidity of surface rock. Assuming this explanation, we are led to attribute to surface 
rocks an average rigidity approximately equal to 6xlO n G.G.S. units, and to the 
Earth as a whole the much higher mean rigidity 1*38 x 10 12 C.G.S. units; further, 
since the ratio of velocities of the first and second set of tremors is approximately 
2 : 1, we are led to assume for X + 2/jl the value 5*53 x 10 12 C.G.S. units, and for v, or 
/x/(X + 2/x), the value £. By analogy to the " tidal effective rigidity " we may introduce 
the phrases " seismic effective rigidity " and " seismic effective modulus of compression " ; 
and the values of these quantities would be 1*38 xlO 12 and 3*69 xlO 12 C.G.S. units 
respectively. When the value of X + 2/x, for the Earth is taken to be 5*53 x 10 12 , the 
corresponding value of s 2 a 2 is 0*625. The results of § 39 appear to warrant the 
conclusion that the moduluses of elasticity of the Earth in its present state are 
sufficiently great to render a spherically symmetrical configuration completely stable. 

41. In obtaining the above values for \+2/x and /x no account is taken of 
gravitation or initial stress, and it is possible that the most appropriate values would 
be a little different from those found above if gravitation and initial stress, to say 
nothing of heterogeneity of density, could be taken into account. For this reason, 
although a complete solution of the problem of wave-propagation in a gravitating 
planet, even when it is regarded as homogeneous, cannot be obtained, the following 
argument may not be without value : — The equations of vibratory motion of a 
gravitating sphere in a state of initial pressure have been obtained in § 3 above. 
From equations (10) and (11) of § 3 we can deduce the equation 

^ = ^±^jv 2 A + 6,9 2 A + ,9 2 r~|, (103) 

and the three equations of the type 



^ = £. JV 2 OT - £/«— -2 —)\ (104) 

where vs xi w y , ts z denote the components of rotation, so that 

9„ -^_^ ctnti\ 

*** x *"""* "^\ ' <"\ > * * * * • * • * * • * * • lJLv/«JI 

In a general way we can see that the terms which contain s 2 in these equations are 
small compared with the remaining terms ; for, if waves of length L are propagated^ 
V 2 A is of the order L~ 2 A, and s 2 A is small in comparison with this in the order s 2 L 2 , 
which is comparable with L 2 /a 2 , since s 2 a 2 is comparable with unity. It would thus 
appear that the velocities of propagation of the waves are not much affected by 
gravitation and initial stress when the wave-length is small compared with the radius 
of the sphere ; and the conclusion would be applicable to superficial waves as well as 
to waves of dilatation and waves of distortion, because such waves are, in any case, 
to be investigated by means of equations of the types (103) and (104). 
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42. In the case of waves of dilatation the argument can be put in a more definite 
shape. Let us suppose that, near a place, the waves a,re plane, so that A is a function 
of x and t P and let us write 

V 1 a = (X + 2/i)//> , . . (106) 

so that V x is the velocity of waves of dilatation when gravitation and initial stress 
are disregarded. We have the equation 

3 2 A , 2 3A , a 2A 1 3 2 A 
ax 2 ox Vi at 

or 

8 



2 > 



dx' 



(e^A) + S 2 (£ - \ s V) ( e ^A) = i- 2 1 (e^A). 



In considering the passage of waves near a place, we may treat the term -~\ix 2 in 
the coefficient of e^' 2 A as a constant ; and then the equation is satisfied by putting 

6^ 2 A = B cos {2rrL~ 1 (x-x Q -Y\t)} > 
provided that 

VV = V/{l-^(^I S v)} (107) 

Since the greatest values of s 2 x 2 are comparable with unity, the value of V\ } the 
local velocity of transmission, is a little less than V 1? or the actual value of \+2jjl is a 
little greater than the seismic effective value. The result (107) may be accepted as 
being not far from the truth in a region large compared with the wave-length, and 
small compared with the radius, and situated at a considerable distance from the 
source of disturbance. 

43. Since the equations of type (104) contain the dilatation as well as the 
components of rotation, it appears that the customary kw of independence of waves 
of dilatation and waves of distortion ceases to hold when gravitation and initial stress 
are taken into account. It appears also that thfe velocities of propagation, both oi 
those waves which are mainly dilatational and of those which are mainly distortional, 
depend on the wave-lengths, and, for the same wave-length, they vary from place to 
place. When the theory can be developed further* these results may possibly prove 
to be useful in explaining the observed irregularities in the propagation of the 
tremors which are recorded in the d&se of great earthquakes. The high values which 
seismic observations lead us to attribute to the elastic constants of the earth as a 
whole are in accord with Lord Kayleigh's view # that great initial stress increases 
the effective values both of resistance to compression and of rigidity. 

* Loc, cit,, ante, p. 173 e 
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PART II. 

A Past State of Gravitational Instability as a Reason for the existing Distribution 

of Land and Water. 

44. Although the conclusion reached by Jeans (1903), that a spherical planet of 
the same size, mass and elasticity as the Earth, in its present state, would be in 
a condition of gravitational stability, is confirmed and strengthened by the present 
investigation, it by no means follows that the Earth has always been in such a state 
as it is now. The fact that the mean density of the Earth as a whole is greater than 
the average density of surface rocks points to a concentration of mass towards the 
centre, and suggests that such a concentration may have come about through the 
elasticity having once been too small for a homogeneous state to be stable. We have 
seen that this would have been the case if the modulus of compression was once as 
small as, or smaller than, that of glass. But we also saw reason to think that, if the 
mean modulus of compression was once decidedly less than that of glass, spherically 
symmetrical states of aggregation would also have been unstable, and the body would 
have existed in some other state. Further, we saw that, if the body was at rest, the 
state in which it would have existed is that which we have described as a state of 
lateral disturbance with a hemispherical distribution of density. The excess of 
density in one hemisphere and defect in the antipodal hemisphere would have existed 
alongside of the concentration of mass towards the centre. 

45. In the paper already cited Jeans (1903) struck out the idea that the 
distribution of land and water on the surface of the globe is associated with a past 
state of gravitational instability. He had found that such instability would manifest 
itself in what has been called above a hemispherical distribution of density. When 
the square of the irregularity is neglected, the figure of a planet at rest, with such 
a distribution of density, is a sphere, but the centre of figure does not coincide with 
the centre of gravity. On taking account of the square of the irregularity, Jeans 
found that the surface of the planet, still supposed to be at rest, would be such as 
can be described roughly as a nearly spherical ellipsoid of revolution, with one half 
slightly flattened at the middle, and the other slightly tapered in the antipodal 
direction. The figure was described as " pear-shaped," the "pear" having a blunt 
end, a sharper end, and a waist. The waters of the ocean would presumably collect 
in the hollow of the waist, and Jeans pointed out that there is some resemblance of 
the shape of the Earth to this figure, although the " stalk" end of the "pear " was 
difficult to discover. 

In the same year a paper was published by W. J. Sollas, # in which it was 
concluded from a discussion of the geographical facts that the shape of the Earth 

* "The Figure of the Earth," < Quart. J. Geol. Soe.,' 59 (1903), p. 180. 
VOL. CCVIL— A. 2 F 
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resembles that of a "pear"; but Sollas' and Jeans' " pears'' have little in common 
beyond the name. Jeans' ideal distribution would consist of a hemisphere which is 
nearly all land, and an antipodal hemisphere which is nearly all ocean, with a central 
island in the middle of this ocean. Sollas' account of the actual distribution is that 
in one hemisphere there is a central continent (Africa) nearly surrounded by a belt of 
seas, while in the antipodal hemisphere there is a central ocean (the Pacific) nearly 
surrounded by a ring of land, the belt and ring being broken at three places, which 
are distributed nearly symmetrically around the centres of the two hemispheres. 
This description suggests very strongly a mathematical account expressed in terms of 
surface harmonics of the third degree. 

If we neglect the rotation of the planet, and regard it as at rest under no external 
forces, we can reach no other result than that reached by Jeans, viz., that, if the 
modulus of compression was once so small that a spherically symmetrical state of 
aggregation would have been unstable, the state of the planet would have been one 
of lateral disturbance with a hemispherical distribution of density. We should not 
be in a position to account at all for the geographical facts as presented by Sollas. 

46. The Earth is a rotating globe, and it is now generally believed to be the 
larger of two fragments into which a single body has been broken up ; the other 
fragment is the Moon. In the early history of the Earth-Moon system the two 
fragments rotated, nearly as a single rigid body ; the period of revolution of the 
Moon was nearly the same as the period of rotation of the Earth. We wish to trace 
the consequences of supposing that the average elasticity of the material was once 
much smaller than it is at present — -that the average modulus of compression was 
more of the order of that of mercury, or even water, than of that of glass or steel, 
and the average rigidity was smaller in comparison with the modulus of compression 
than it is to-day. We have the problem of determining the distribution of density 
within the planet, and the consequent shape of its surface. The problem cannot be 
solved completely, but we can make some progress with it ; and we can then attempt 
to discover the extent to which our results accord with geographical observation. In 
so far as the accord is good we may regard geography as supporting the hypothesis 
as to the past state of the Earth. 

Illustration of the Nature of a Hemispherical Distribution of Density. 

47. We have reason to think that, in the absence of rotation and external forces, 
the planet, if of sufficiently small elasticity, would have been in the state which we 
have described as a state of lateral disturbance with a hemispherical distribution of 
density. Before proceeding to take account of the rotation and external attractions, 
we consider further the nature of such a disturbance. For this purpose we take the 
problem of a spherical body, homogeneous when unstrained, and devoid of all rigidity, 
and suppose that in the initial state the s§lf- attraction of the body is balanced by 
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hydrostatic pressure. We suppose also that the law of elasticity of the body is that 
the increment of pressure is proportional to the increment of density. We show 
that equilibrium is possible in strained states, in which the excess of density at any 
assigned distance from the centre is proportional to a spherical surface harmonic of 
the first degree. 

In the initial state the pressure p Q and potential V are given by the formulse 

Po = inypo 2 (a 2 -r 2 ), V = %iryp (3a 2 -r 2 ). 

In the strained state the pressure P, density p, and potential V are expressed by 

the formulae 

P=p +X£ a> = />o(1 + £), V = V +W, 

where £ denotes the condensation. The equations of equilibrium are 

dx dx ty ty ' dz dz J 

and W is connected with £ by the equation 

V 2 W = - iirypol 

When terms of the second order in the small quantity ^ are neglected, the 
equations of equilibrium become three equations of the type 

^ _|7rw ° 2 ^~ x i =0; • (108) 

and, on eliminating W, and writing s 2 for $iryp 2 j\ we have 

V*£+£r^ + W£=0 (109) 

This equation is satisfied by putting 

where A is an arbitrary constant and o) x is a spherical solid harmonic of the first 
degree, and this is the most general form of solution in which £ is finite at r = 0, and 
is proportional when r = const, to a surface harmonic of the first degree. The 
additional potential W has the form 

W = 47ryp {lA^V^ 2 o) 1 + F 1 }, 

where ¥ x denotes a spherical solid harmonic of the first degree. 
Let the bounding surface become 

r = a + U a . 
2 F 2 
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Since the pressure vanishes at this surface, the expression 



2. 
3 



7ry P o 2 {a 2 ~(a + U,) 2 } + XA ( 1 -~}s 2 a 2 ) e~*** a ( Wi / r ) 



vanishes, or we have, neglecting U/ ? 

so that TJ a contains the same surface harmonic as ojj. The form of ¥ x is determined 
by the condition that W is the potential of a distribution of density p £ through the 
volume of the sphere r = a, together with a distribution of density p JJ a on its 
surface. Just as in § 14, this condition leads to the equation 
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48. Now let the bounding surface in the strained state be 

r = a+b cos 6, 

which represents a sphere with its centre at a small distance b from the origin in the 
direction of the axis of the harmonic. We find 



f = 
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If s 2 a 2 > 5, the condensation is greatest near the centre, and it is positive on the 
side remote from that towards which the surface is displaced, so that the centre of 
gravity is displaced in the opposite sense to the surface. The distance of the centre 
of gravity from the origin is easily proved to be 5b/(s 2 a 2 —5). 




Axis of Harmonic 



Fig. 1. 




Fig, 2. 
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The variation of the excess density along the axis of the harmonic is illustrated in 
fig. 1. The surface r = a (1 + € cos 9) can be an equipotential surface if 

ae = — bbj(s 2 a 2 — 5), 

and thus a sphere of radius a with its centre at the displaced centre of gravity is an 
equipotential surface. The relative situation of the bounding surface and of this 
equipotential is illustrated in fig. 2, in which denotes the undisplaced centre, C the 
centre of the displaced surface, and G the centre of gravity of the strained sphere. 
The figures are drawn for the case in which s 2 a 2 =10. 

The type of disturbance which has been called above a lateral disturbance with a 
hemispherical distribution of density would be the same in a body possessing some 
degree of rigidity, but the numerical details would be different. 

49. If the equipotential surfaces of a nearly spherical body, with a nearly 
symmetrical distribution of density, are referred to the centre of gravity of the body 
as origin, their equations take such forms as 

r = a(l + € 2 S 2 +€ 3 S3+...), 

in which e 2 , ... denote small coefficients, and S 2 , ... denote spherical surface harmonics 
of degrees indicated by the suffixes. There is no term of the form efi^ In the case 
of the Earth, the coefficients € 2 , ... can be determined by means of pendulum 
experiments. If we referred to a different origin, near the centre of gravity, a term 
of the form € 1 S 1 would be introduced, but the coefficient e ± could not be determined by 
means of pendulum experiments, for it does not affect the formula for the variation 
of gravity over the surface. # If we choose an origin in accordance with geometrical 
considerations, e.g., as the centre of that oblate spheroid which most nearly coincides 
with the surface of the ocean, the results of pendulum experiments cannot tell us 
whether this origin coincides with the centre of gravity or not. 



Effect of Rotation upon a Planet with a Hemispherical Distribution of Density, 

50. In all the preceding work the rotation of the Earth has been neglected. We 
have now to consider the effect of rotation upon a nearly spherical planet which, in 
the absence of rotation, would have a hemispherical distribution of density. To 
simplify the analysis, we shall disregard the concentration of mass towards the 
centre and also the rigidity of the body. We shall take as the " initial " state of the 
body a state in which the density is uniform and the stress is hydrostatic pressure, 

* The result may be inferred from Stores' investigation of the " Variation of Gravity over the 
Surface of the Earth," Cambridge, ' Trans. Phil. Soc.,' 8 (1849), or ' Math, and Phys. Papers,' vol. 2, 
Cambridge, 1883. It is easy to prove it independently. 
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while the body is rotating, as if rigid, about the axis of z with angular velocity co ; 
and we shall seek a strained state in which the body could exist without the 
application of any external forces, this state being such that, in the absence of 
rotation, the distribution of density would be hemispherical. In the notation of § 47 
the equations of steady motion of the body are 



2 av ap 

pcox = p- — , 

ox ox 



2, 



posy 



av^ap 

dy dy ' 



av 3P 



= p^-^. . . (110) 

0% oz 



The initial state is determined by the same equations with p^ V , p Q substituted for 
p, V, P. Now the initial figure is an oblate spheroid, and the initial form of V is 

V = const. -J {A' (x 2 +^ 2 ) + CV} ? . 

where A' and C are constants ; also the initial form of P is 

p Q = const. + p { V + ^co 2 {x 2 +y 2 )}, 

= const. - ^p { (A' - co 2 ) (x 2 + y 2 ) + C f z 2 } . 

When we write, as in § 47, 

and neglect terms which cancel on account of the values of p and V , and 
also neglect terms which are of the second order in the small quantity f, the 
equations (110) become 



p Q o) 2 xt; = -pgA/x + p -z— -X^. , 

ox ox 

' PQO)yg = ~ M A 'y + /On -r— ~ a 



*\ 



^ ? 



dy dy 

= — p Q gC f Z + pa-^~~-\ b 



^ • • • • • • ( 1 i I ) 



3z 



a» • 



y 



Now we have the equations 



2A' + C - 47r r/ ) , V 2 W = -47ry/> £ 
and therefore we can eliminate W and obtain the equation 

vf + (w-we)f + ft{(A'-rf)(«| + »«) + (y,|} = o, . (ii 2 ) 

where s 2 is written for ^iTyp 2 /k. If to were zero, A/ and C would both be equal to 
f 7ry/) , and we therefore put 



A' - iirypo + A", C - iTtypo + C" ; 



GRAVITATIONAL STABILITY OP THE EAETH. 223 

then equation (112) becomes 

V^ +s V^+6^=2a, 3 ^+^{(A"- w 2 )^||+ 2 /^)+C // 2||}. . (113) 

The left-hand member of this equation is the same as that of equation (107) in § 47 
above ; and therefore, when co 2 is neglected, f can be of the form £ l9 where 

£ = A(l-i S V) e -^V, 

the notation being the same as in § 48. Now we shall suppose that <o 2 is not large, 
so that we may treat f x as an approximation to £ and substitute & for £ in the 
right-hand member of equation (113), for all the terms of this member are small of 
the order o> 2 £ We are then neglecting £ 2 , but not <o 2 £ To obtain a second 
approximation, we put 

£ = ii + f ', 

and seek a particular integral of the equation 

V 2 g' + s 2 r^-+6s 2 e = 2a> a £ £+£ j(A"- w 2 )(^ +yM] +C^^1. (114) 

There would be no special difficulty in obtaining a solution of the equation, but it 
will be sufficient for our purpose to find the form of the solution. The function £ 
may be expressed in terms of polar co-ordinates r, 0, $ in the form 

^ =f(r) r (<% sin cos <£ + & sin # sin (fy + yj cos 0), 

where a u /3 U y x are constants, and f(r) is a certain function of r which has been 
determined. Hence we have 

x~ +i/~ = rfMfo^ sin cos <& + & sin sin 6) 

dx dy ; 

+ r 2 f (r) sin 2 (a 2 sin cos </>-f j8 x sin sin <£+yi cos 0), 
% Si = rf(r) y x cos + r 2 f r (r) cos 2 6 (a x sin cos <J> + j3 x sin sin <£ + yi cos 0) ; 

C//V 

and these can be expressed in the forms 

x k^ +Vjr- = r {/( r ) +i r f f ( r )} sin-0 (a x cos <£ + /3i sin <£) +fr 2 / / (r) y x cos 

—r 2 f^(r)(co$ 2 0—j-) sin (a x cos <£ + /3i sin <£) —r 2 f f (r) y x (cos 3 0— f cos 0), 

zSl z=z^r 2 f'(r) sin (a x cos <£+/3i sin <jf>) + ^{/(^) +f'^ / (' r )}yi cos 

+ ry v (r) sin (cos 2 0— §-)(«i cos <£ + /?i sin <£) + r 2 f , (r) y x (cos 3 0— f cos 0) 

Hence the right-hand member of (114) can be expressed as a sum of terms each of 
which is the product of a function of r and a spherical surface harmonic, ajxd the 
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surface harmonics which occur are those of the first degree and the following 
harmonics of the third degree : — 

cos 3 0— f cos #, (cos 2 9— \) sin 6 cos <f> 9 (cos 2 0—j) sin 6 sin <f>. 

To each of these terms there corresponds a term of the same form in £', and 
therefore also in P, or p + ^; and it follows that the displacement of the bounding 
surface from its initial form (which is a slightly elliptic oblate spheroid appropriate to 
the rotation) is expressed by a radial displacement;, which consists of a part propor- 
tional to a spherical surface harmonic of the first degree, together with parts propor- 
tional to the above surface harmonics of the third degree. In like manner all the 
terms of the additional potential W are the products of functions of r and surface 
harmonics, which are either of the first degree or are the above harmonics of the 
third degree ; but the coefficients of the various harmonics in W are different from 
their coefficients in £ The equation of the equipotentials under gravity, modified by 
the rotation, is 

V e + W + j^co 2 (x 2 + y 2 ) = const. , or p /p + W = const. ; 

and thus the situation of the bounding surface relative to the equipotentials is 
expressed by a difference of radii at corresponding points, this difference being a sum 
of terms of the form 6S, where b denotes a constant and S denotes a surface 
harmonic ; and the surface harmonics which can occur are those of the first degree 
and the three of the third degree written above. 

51. It appears from this investigation that, if a gravitating body, which is rotating 
about an axis, has so small a modulus of compression that, if the body were at rest, a 
spherically symmetrical distribution of density would be unstable, it would tend to take 
up a state in which the distribution of density would not be exactly hemispherical, 
but the excess density would also contain terms expressed by spherical harmonics of 
the third degree. The figure of the body would differ from the oblate spheroidal 
figure appropriate to the rotation by a radial displacement at each point ; and this 
displacement would be expressed partly by spherical surface harmonics of the first 
degree, indicating that the centre of gravity does not coincide with the centre of 
figure, and partly by spherical harmonics oi the third degree. If the body were 
entirely devoid of rigidity, the oblate figure appropriate to the rotation would be the 
same as that of an equipotential surface under gravity, modified by the rotation ; and 
the figure of the body, as determined by difference of level above or below a certain 
equipotential surface, would be an harmonic spheroid of the third degree, and the 
situation of the body would be that of such a spheroid when displaced towards one 
side. If the body possessed some rigidity, the oblate figure appropriate to the 
rotation would differ a little from that of a nearly coincident equipotential surface, 
and the shape of it, determined as before, would be that derived from a certain oblate 
spheroid of small ellipticity by a displacement proportional to a surface harmonic of 
the third degree. The surface harmonic would be of a somewhat specialised type. 
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Effect of certain External Forces. 

52. The effect of forces such as the attraction of the Moon at the time when its 
period of orbital revolution did not differ much from the period of rotation of the 
Earth would be to draw the planet out into a shape more nearly ellipsoidal, with 
three unequal axes, than spheroidal. If the planet could have had a symmetrical 
shape it would have been practically ellipsoidal, and the surfaces of equal density 
would have been ellipsoids. Whereas the effect of rotation is the same as that of 
forces derived from a potential of the second degree, symmetrical about the axis of 
rotation ; such forces as we are now considering are derived from a potential, the 
most important terms of which would be of the second degree, but not symmetrical 
about the axis of rotation. If the elasticity was too small for an ellipsoidal figure to 
be stable, the planet would have been in a disturbed state, the nature of which can 
be inferred from the preceding investigation. We have only to replace in § 50 the 
initial potential, modified by the rotation, by a general expression of the second degree 
in the co-ordinates. The only change that would be made in the result would be that 
those terms in the radial displacement which are expressed by harmonics of the third 
degree would not be of the specialised type introduced by the rotation, but would be 
of general type. The figure of the planet would be derived from the ellipsoidal 
figure appropriate to the rotation, and to the external forces, by a radial inequality 
expressed by surface harmonics of the first and third degrees. The equipotential 
surfaces would be obtained from the ellipsoidal equipotentials appropriate to gravity, 
modified by the rotation and the external forces, by surface harmonics of the same 
degrees. The result would be that the shape of the planet, as determined by 
difference of level above or below a certain equipotential, would be a wrinkled 
ellipsoid, displaced towards one side ; and the wrinkle would be expressible by means 
of a spherical surface harmonic of the third degree. 

The Problem of the Shape of the Lithosphere. 

53. The problem of determining the form of the equipotentials near the surface of 
the Earth includes the problem of determining the figure of the surface of the 
ocean (the "hydrosphere"). The equipotentials which lie outside the nucleus (or 
" lithosphere") on one side, and sufficiently near to it, cut the surface of the 
lithosphere towards the other side. Among these equipotential surfaces that one 
which, outside the lithosphere, coincides with the surface of the ocean is known as 
the " geoid." The surface of that part of the lithosphere which lies outside the geoid 
is occupied by land, and can be observed directly ; the surface of that part which lies 
within the geoid can only be observed indirectly by means of soundings. We have 
no means of investigating the form of the surface of this part of the lithosphere 
except by estimating its depth at a point below the geoid. The most important 
deviations from sphericity both of the lithosphere and of the geoid are of such 
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a nature that these surfaces are nearly oblate spheroids. If the lithosphere were 
exactly in the form of an oblate spheroid, and its centre of gravity coincided with its 
centre of figure, it would either lie entirely within the geoid or would protrude from 
it symmetrically at the North and South Poles. Owing to the rigidity of the 
lithosphere the ellipticity produced in the geoid by rotation would be slightly greater 
than that produced in the lithosphere, and thus there is a tendency to lay bare the 
polar regions ; but, since the land of the globe does not consist of two circular islands 
at the poles, there are other deviations from sphericity, both of lithosphere and geoid, 
and the relative amounts of these at different places can be expressed by the difference 
of radii drawn from the centre of gravity. According to the theory which has been 
here advanced this difference of radii should be, at least in its general features, 
expressible as a sum of spherical harmonics of the first, second and third degrees. 

54. It is easy to verify the presence of some of these harmonics. The effect of 
a term of the first degree would be to make the lithosphere protrude from the geoid 
towards one side. If this term were the only one, the land of the globe would form 
a circular island or continent. It is the fact that most of the land is in one 
hemisphere. The great circle of the globe which contains most land has a pole 
situated between Orleans and Le Mans # (latitude 48° N., longitude 30 ; E.). Again, 
the zonal harmonic of the third degree vanishes at three circles, one being a great 
circle. If this term were the only one, the land of the globe would consist of 
a circular island surrounded by a belt of ocean in one hemisphere, and in the 
antipodal hemisphere there would be a circular ocean surrounded by a ring of land. 
This arrangement corresponds to two features of Sollas' description of the Earth's 
surface. The nearly symmetrical breaking at three places of the belt and three of 
the ring, which he also noticed, indicates the presence of the sectorial harmonic of the 
third degree. If we refer to the polar axis, instead of any other morphological axis, 
the presence of the zonal harmonic of the third degree is indicated by the existence 
of an Antarctic continent, and by the fact that most of the land of the globe is north 
of the Equator. The harmonic of the third degree and second rank, referred to the 
polar axis, vanishes at the Equator and at four meridians symmetrically placed. If 
this term were the only one, then, in two northern quadrants there would be land, 
and also in the two alternate southern quadrants, an arrangement which suggests 
Central Asia and North America as the land quadrants of the northern hemisphere, 
Australia and South America as those of the southern. 

Spherical Harmonic Analysis of the Distribution of Land and Water, 

55. By such arguments as the foregoing, and by some trials with small numerical 
coefficients for the various harmonics, I had convinced myself that many features of 
the distribution of land and water could be represented by means of harmonics of the 
third degree, when Professor H. H. Turnee suggested to me the advisability of 

* E, ButiOKNER, 'Die feste Erdrinde und ihre Formen/ Wien, 1897. 
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adopting a systematic process for the discovery of appropriate coefficients. He very 
kindly made, and placed at my disposal, a rough preliminary calculation, and the 
results were sufficiently encouraging to warrant the undertaking of a considerable 
piece of computation. A professional computer was employed for a time, but 
eventually I relied upon my own calculations, taking many precautions to ensure 
accuracy. The systematic process consists in devising a function to represent the 
" value of land" at any point, and determining, by the method of approximate 
quadrature, the coefficients of an expansion of the function in spherical harmonics. 
The results of such a computation clearly depend upon the chosen " value of land," 
and judgment must be exercised in selecting appropriate values. Little importance 
can be attached to the heights of mountains, because the highest mountain ranges 
are, geologically speaking, modern, the ancient mountains being worn down by 
denudation and erosion. Too much importance is not to be attributed to the actual 
coast-line, because this line is subject to many causes of change. The coast-line is 
but one of the contour-lines of the continental block (the geoid being the level of 
reckoning), and the shape of the block at considerable depths differs a good deal from 
that at the surface. At mean-sphere-level (8400 feet below sea-level) the continents, 
with the exception of the Antarctic continent, form, a continuous block. # The Arctic 
Ocean is reduced, so far as is known, to a trough running nearly along the meridian 
of Greenwich, from about latitude 65° N. to about latitude 80° N. It may extend to 
the North Pole and surround it. The polar block spreads southwards in two great 
masses — America and Eurasia. These are joined through the British Isles, Iceland 
and Greenland on the one side, and across Behring's Strait on the other ; the contour- 
line at mean- sphere-level runs practically along the 60th parallel between America 
and Europe and along the 50th parallel between America and Asia. The Eurasian 
division of the block forks near the Persian Gulf, and tapers southwards in two 
branches, one containing Africa and the other the Malay Peninsula, adjacent islands, 
Australia and New Zealand. The Red Sea does not go down to mean-sphere-level, and 
the Mediterranean does so only in two small patches. The American division of the 
block is continuous across the Gulf of Mexico, the West Indies and the Caribbean 
Sea, which, at this depth, equally with Mexico, Central America, and the Isthmus of 
Panama, form part of the ridge joining North and South America. The ridge has 
some local depressions. The block tapers towards Cape Horn, in the neighbourhood 
of which, however, it has a great eastward extension, and this extension turns 
westward and nearly joins the northern continental block to the Antarctic continental 
block through the South Shetland Islands.*!" The Antarctic block also shows a 

* The information here detailed in regard to the distribution of the continental blocks and oceanic 
regions at mean-sphere-level is taken from a map drawn by H. R. Mill in ' The Scottish Geographical 
Magazine ' . (Edinburgh), vol. 6 (1890), p. 184. Reference may be made to the rough map on p. 237 below. 

t It is now known that the depth of the channel is not so great as it was for a long time supposed to be. 
See a paper by W. S. Bruce in 'The Scottish Geographical Magazine' (Edinburgh), vol. 21 (1905), p. 402. 
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northward extension towards Australasia. The contour-line of the continental blocks 
at mean-sphere-level is a very important and fairly well ascertained datum of the 
problem. If, however, we attend exclusively to it, we are liable to emphasise unduly 
those parts of the block which do not rise above the level of the sea. 

56. I calculated the coefficients of a spherical harmonic expansion up to harmonics 
of the third degree for two different assumptions as to the " value of land." . In the 
first assumption the value — 1 was attached to those points of the surface which are 
below mean- sphere-level and the value to those points which are above it. In 
the second assumption the value 1 was attached to those points of the surface which 
are above sea-level and the value to those below it. The coefficients obtained by 
the two assumptions were then added. The somewhat greater importance of the 
mean sphere may perhaps be sufficiently represented by the result that the maxima 
obtained by using the first set of coefficients are larger than those obtained by using 
the second set. The combined distribution for the two sets of coefficients is shown in 
the following table, in which 6 stands for co-latitude measured from the North Pole, 
and <f> for longitude measured eastwards from the meridian of Greenwich : — 
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1 


1 


1 


-< 

i 


1 




-i 
- i 








125 


-1 


-1 


-1 


-1 




1 


1 










1 


1 




-1 


-1 






130 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 






i 


-1 


!-l 


-1 




1 


135 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 






-1 


-1 


-1 








140 


1 


~1 


-1 


- 1 


-1 


-1 


-1 


-1 


-1 




-1 


-1 


-1 


-1 


i-i 


-1 


-1 


-1 


145 


„1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-l 


-1 


-1 


-1 


150 


-1 


1-1 


-1 


-1 


-1 


-1 


-1 


-1 


— i 


-1 


_ i 


-1 


-1 


-1 


-l 


-1 


-1 


-1 


155 






1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


! i 


1 


i 








160 






1 


1 


1 


i 
i 


1 


1 


1 


1 


1 


1 


i 


1 


l 








165 


1 


1 


1 


1 


i 


i 


1 


1 


X 


1 


1 


1 


i 


1 




1 


j 




170 


1 


1 


1 


1 


1 


i 


1 


X 


1 


1 


1 


1 


i 


1 


i 


1 


1 


1 


175 


1 


1 

i 


i x 

1 


1 


1 


i 


1 


1 


1 


1 


1 


1 


i 


1 


! i 

i 


1 


1 


1 

i 
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Table (continued). 



</> 



180°. 



185°. 







5 

10 

15 

20 

25 

30 

35 

40 

45 

50 

55 

60 

65 

70 

75 

80 

85 

90 

95 

100 

105 

110 

115 

120 

125 

130 

135 

140 

145 

150 

155 

160 

165 

170 

175 



-1 



1 
1 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



1 
1 
1 
1 



1 
1 



-1 



1 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



1 
1 



190°. 


195°. 


-1 


-1 




1 


1 


1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


i 
— ± 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


_1 


-1 


-1 


-1 


-1 


1 


1 


1 


1 

i 



200°. 



-1 



1 
1 
1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



1 
1 



205°. 



-1 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



1 
1 



210°. 



-1 



1 
1 
I 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 



1 
1 



215°. 



-1 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



1 
1 



220°. 



1 
1 
1 

-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 



1 
1 



225°. 



-1 -1 



1 
1 
1 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



1 
1 



230°. 


235°. 


240°. 


-1 


-1 


-1 
1 


1 


1 


1 


1 


1 


I 


1 


1 


1 


1 


1 


1 
1 
1 
1 


-1 






-1 


-1 




-1 


-1 


-1 


-1 


-1 


-1 


--1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


- 1 


-r 


-1 


-1 


-i 


-1 


- 1 


-l 


-1 


-1 


-l 


-1 


-1 


-l 


-1 


-1 


-i 


-1 


-1 


-i 


-1 


-1 


-i 


-1 


-1 


-l 


-1 


-1 


-l 


-1 


-1 


-l 


-1 


-1 


-l 


-1 


1 


l 


1 


I 


l 


1 



245° J 250°. 



1 

1 
1 
1 
1 
1 
1 
1 
1 
1 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



1 
1 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



1 
1 



255°. 



-1 

1 
1 
1 
1 
1 
1 
1 
1 
1 
1 



-I 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 
-1 



1 
1 



260°. 


265°. 


-1 


-1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 






1 


-1 


-1 


-1 


-1 


-1 


- 1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


... 1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


1 


1 


1 


1 

i 
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Table (continued). 



\ * 






































\ 


270°. 


275°. 


280°. 


285°. 


290°. 


295°. 


300°. 


305°. 


310°. 


315°. 


320°. 


325°. 


330°. 


335°. 


340°. 


345°. 


350°. 


355°. 


o \ 






































o 

5 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


10 








1 










1 


1 


1 


1 


1 












15 


1 


1 


1 




1 


1 


1 


1 


1 


1 


1 


1 


1 












20 


1 


1 


1 


1 












•1 


1 


1 


1 












25 


1 


1 


1 


1 












1 


















30 








1 










-1 








-1 


-1 










35 








1 


1 


1 






-1 


-1 


-1 


-1 


-1 


-1 


-1 








40 


1 


1 


1 


1 


1 


1 


1 






-1 


-1 


-1 


-I 


-1 


-1 






1 


45 


1 


1 


1 


1 


1 


1 




-1 


-1 


-1 


-1 


-1 


-1 


-I 


-1 


-1 




1 


50 


1 


•1 


1 


1 




-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


1 


1 


55 


1 


1 


1 




-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 




1 


60 


1 








-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 




1 


1 


65 












-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 






1 


1 


70 














-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 




1 


1 


1 


75 


1 












-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 




1 


1 


1 


80 




1 






1 






-1 


-1 


-1 


-1 


„1 


-1 


-1 








1 


85 


-1 


-1 




1 


1 








-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 






90 


-1 


-1 




1 


1 


1 


1 


1 


1 




-1 


-1 


-1 


-1 


-1 


-1 


-1 


-i 


95 


-1 




1 


1 


1 


1 


1 


1 


1 


1 


1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


100 


-1 


-1 




1 


1 


1 


1 


1 


1 


1 


1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


105 


-1 


-1 






1 


1 


1 


1 


1 


1 




-1 


-1 


-1 


-1 


-1 


-1 


-1 


110 


-1 


_ i 

jl. 


-1 




1 


1 


1 


1 


1 


1 




-1 


-1 


-1 


-1 


-1 


-1 


-1 


115 


-1 


-1 


-1 




1 


1 


1 


1 


1 




-1 


- 1 


-1 


-1 


-1 


-1 


-1 


-1 


120 


-1 


-1 


-1 




1 


1 


1 


1 


1 




-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


125 


-1 


-1 


-1 




1 


1 


1 


1 




-1 


-1 


-1 


-1 


-1 


-1 


-I 


-1 


-1 


130 


-1 


-1 


-1 




1 


1 








-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


135 


-1 


-I 






1 


1 




-1 


-1 


-I 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


140 


-1 


-1 






1 






-1 


-I 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


145 


-1 


-1 






















-1 


-1 


-1 


-1 


-1 


-1 


150 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


-1 


155 


















-1 


-1 


-1 


-1 


-1 


-1 


-I 


-1 


-1 


-1 


160 


1 


1 


1 


1 






















-1 


-1 


-1 


-1 


165 


1 


1 


i 


1 


1 


1 


1 


1 


1 


1 


















170 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 










175 


1 


1 


1 

1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 


1 
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57. The surface harmonics of the first degree expressed in ordinary spherical polar 

co-ordinates 0, 6 are • n , • a • ± a 

9 ^ sin v cos 9, sin u sin <p, cos u ; 

and any spherical surface harmonic of the first degree can be expressed in the form 

(p cos <J> + q sin <£) sin 6-\-r cos 0, (115) 

where p, q, r are numbers. The spherical surface harmonics of the second degree 

sin 20 cos (f>, sin 20 sin <f> 9 sin 2 cos 2(f), sin 2 sin 2(f), 3 cos 20+ 1 ; 
and any spherical surface harmonic of the second degree can be expressed in the form 

(a cos <£ + /3sin <£) sin 20+ (y cos 2(£ + S sin 2<£) sin 2 + e(3 cos 20+1). . (116) 

The spherical surface harmonics of the third degree are 
(i.) The zonal harmonic fi CQg3 ^ 3 cos ^ . 

(ii.) The tesseral harmonics of the first rank 

(5 cos 2 0—1) sin cos <f> 9 (5 cos 2 0—1) sin sin <f> ; 

(iii.) The tesseral harmonics of the second rank 

sin 2 cos cos 2</>, sin 2 cos sin 2(f> ; 

(iv.) The sectorial harmonics 

sin 3 cos 3<jf>, sin 3 sin 3<£. 
Since 

5 cos 3 0— 3 cos = f (cos 30+f cos 0), 

(5 cos 2 0— 1) sin = i (sin 0+ 5 sin 30), 

sin 2 cos = -J (cos 0— cos 30), 

sin 3 = -J (3 sin 0-- sin 30), 

any spherical surface harmonic of the third degree can be expressed in the form 

aW + X (b cos (f> + c sin <£) + Y(cl cos 2</> + e sin 2</>) + Z(/cos 3(f)+g sin 3<£), . (117) 

where a, &, c, d, e, f 9 g are numbers, and 

W = cos 30+ (0*6) cos 0, 
X = sin 0+5 sin 30, I 

Y = cos 0- cos 30, ^ (118) 

Z = 3 sin 0— sin 30. 



j 



* The form 3 cos 20+ 1 for the zonal harmonic is 4 (f cos 2 #- 1~) 5 and is taken as being more convenient 
for calculation. 
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58. Let F(0, <f>) denote the function to be expanded. The coefficients are expressed 
by equations of the type 



•2tt 



'IT 



•2rr 



p d<f> d0(cos <j> sin 0) 2 sin = dj>\ dffF(0 9 <£)(cos <f> sin 0) sin 0. . (119) 



Jo 



The factors multiplying the coefficients p, &c. in the left-hand members are the 
integrals of the squares of the several harmonics over the surface of a unit sphere. 
The integrals in the right-hand members are the integrals, over the surface of the 
same sphere, of the product of the function to be expanded and the corresponding 
harmonics. The values of the integrated squares multiplying p, &c., are recorded in 
the following table : — 



Coefficient. 


Value of integrated 
square. 


Reciprocals. 


P, £ r 


iir 
<) 

■UtT 

.1 o 


( W) x -AV^r- 1 

\ o o / 2 t> o 


a 3 /J, y, 8 


(¥-) x AV^ 1 


e 


(40 x AV^- 1 


a 


1 t o 


(5) x ¥ ^- J 


b, c 


5 12_ 
21 " 


(0 • 3) x Tp^-i 


d, e 


25.6^ 
10 5 /r 


(3) x ¥ 3 JVx-i 


f> 9 




(0-5)x^-i 



Since the ratios only are relevant, the integrals in the right-hand members of such 
equations as (119) are to be multiplied by the numbers in brackets in the third 

column. 

59. To evaluate integrals of the type in the right-hand member of (119), when the 
value of F (#,<£) is given by the table of §56, or any similar table, we treat the 
integral as a double sum, e.g., 



m=7l / _ \ w=35 , _ 



m 



-vx/nTT mir\ ( mrr • mr\ • nir 
- , , _ , , r — , -—- cos - — ■ sin — sin — - : 
=o V36/ »=i \36/ \36 ? 36/V 36 36/ 36 



then we have to evaluate such a double sum as 



S SI ^— cos — r- sm 



m = W = l 



V36' 36/ 



36 



36 
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We sum first with respect to m ; but in forming the sum we take account of the 
fact that sin 2 (^7r/36) does not change when n is replaced by 36— n. For example, 
let F be equal to 1 at the points indicated in the table, and zero at other points. 
Then the contribution to the terms containing any m of the two parallels given by n 
and 36— n is either 0, 1, or 2, according as a 1 occurs on neither parallel (for the 
particular m in question), on one, or on both. This number 0, 1, or 2 is to be 
multiplied by the value of cos(m7r/36) for the chosen m ; but the same value for the 
cosine occurs at the meridian given by 72— m, and the same numerical value with the 
opposite sign occurs at the meridians given by 36 — m and 36 + m. We condense into 
one term the contributions of the eight points given by n, 36 — n, m, 72 — m, 36 ± m > 
and take the ranges of m and n to be respectively to 17 and 1 to 18. Thus, as the 
multiplier of cos (mn/SG) sin 2 (mr/36), we have an integral number which necessarily 
lies between — 4 and 4, and may be zero, and we have transformed the sum. into a 
double sum of the form 

2 2 Jd (n, m) cos — sm 



Wi = 71: 



a v ' J 36 36 



where F' is the number in question. The most troublesome part of the process is the 
determination of F'. When F' has been found it is very easy to form the sum of 
such a series as that written immediately above by summing first with respect to m 
and then with respect to n. When we are dealing with tesseral harmonics of the 
second rank, we can thus condense into one term the contributions of 16 points of 
the table, and, when the tesseral harmonic is of the third rank, those of 24 points. 
Much labour is saved by going through this process, troublesome though it is, and 
much greater accuracy can be secured, because in the multiplication of cos (m7r/36) 
by F', when F' is, say, 5 or 6, and the value of the cosine to any chosen number 
of decimal places is used, it is easier to correct the figure in the last place than it is 
when the same cosine occurs five or six times in a long column of figures which have 
to be added together. 

60. By the use of this method I computed the values of the coefficients p, &c, for 
the function F (#, cj>) which is given by the — 1's in the table of § 56, the 1's being 
replaced by zeros. Up to the stage of summation with respect to m, inclusive, I kept 
four decimal figures. Of the terms of the type 



17 ' mn \ • 9 nrr 



2 F ; (n, m) cos — — ) . sin 



m = 



36/ 36 



I then kept two decimal figures, formed the sums with respect to n, and multiplied 
them by the corresponding numbers placed in brackets in the third column of the 
table in § 58. This process gave the coefficients in the second column of the annexed 
table. The integral parts only were retained. I computed the values of the 
coefficients p, &c, in the same way for the function given by the 1's in the table of 
§ 56, the — 1's being replaced by zero. This process gave the coefficients in the third 

2 H 2 
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column of the annexed table. It is to be understood that in both cases the common 

factors - — -ir"" 1 and (77-/3 6 ) 2 have been omitted, 
256 





F (<9, <£) == or - 1. 


F (0, <£) = 1 or 0. 


Sum. 


J? 


604 


557 


1161 


2 


495 


329 
630 


824 


r 

1 


777 


1407 


a 


350 


243 


593 


P 


295 


366 


661 


7 


-443 


-223 


- 666 


8 


-291 


68 


-223 


€ 


185 


98 


283 


a 


-213 


-134 


-347 


b 


- 73 


- 71 


- 144 


G 


- 29 


38 


9 


d 

e 


-338 


-256 


-594 


396 


351 


747 


f 


56 
203 


26 


82 


9 


122 


325 



The Continental Blocks and Oceanic Regions as expressed by Spherical Harmonics 

of the First, Second and Third Degrees. 

61. I then computed the values of the harmonics expressed by (115), (116), (117), 
for values of and <£, which are multiples of 15° (or y^tt), using first the coefficients 
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given in the second column of the table in § 60, and then the coefficients which are 
given in the third column of the same table. Finally I added the values belonging 
to the same 6 and <£. The results are shown in the diagram (fig. 3), where the fine 



165° 150° 135° 120° 105° 90° 75° 60° 45° 30° 15° 0° 15° 30° 45° 60° 75° 90° I0S° 120° 135° 150° 165° 




!65° 150° 135° 120° !Q5° 90° 75° 60° 45° 30° 15° 0° 15° 30° 45° 60° 75° 90° 105° (20° 135° 150° 165° 



Fig. 3. 



165° ISO 135° 120° 105° 90° 75° 60° 45° 30° !S° 0° 15° 30° 45° 60° 75° 90° 105° 120° 135° 150° 165° 




165° 150° 135° 120° 105° 90° 75° 60° 45° 30° IS 0° 15° 30° 45° 60° 75° 90° 105° 120° 135° 150° 165° 

Fig. 4. 

continuous line is the contour-line along which the calculated harmonic inequality 
vanishes, the heavy continuous line is the contour-line along which this inequality is 
10 per cent, of its maximum below zero, and the dotted line is the contour-line along 



238 PROFESSOR A. E. II. LOYE ON THE 

which this inequality is 10 per cent, of its maximum above zero. It is to be observed 
that an inequality expressed by harmonics of uneven degrees has numerically equal 
values with opposite signs at antipodal points, and therefore the area on the sphere 
within which such an inequality is positive is equal to the area within which it is 
negative. But this equality of positive and negative areas does not hold when the 
harmonics of the second degree are present. A rough calculation showed that the 
zero line of the inequality illustrated in fig. 3 divides the surface of the sphere into 
two unequal areas, and the inequality is negative in the larger area. The excess of 
the negative area above the positive is nearly 10 per cent, of the whole surface. 
The heavy line in fig. 3 corresponds more nearly than the other lines to the 
principle by which geographers construct the contour-line at mean-sphere-level. The 
diagram in fig. 3 suggests many features of the outline of the continental block, and 
there can be no doubt that the coefficients could be adjusted so as to secure a better 
agreement. # It seems best, however, to record the results as they are. For the sake 
of comparison a rough map of the world is added (fig. 4). The heavy continuous line 
is the outline of the continental block at mean-sphere-level, and the fine continuous 
line is the coast-line. I have not attempted to draw the map with minute accuracy ? 
and have omitted many small islands and some small enclosed patches of deep sea, 
because the object aimed at is a comparison of the general features of the map of 
the world with those of the diagram in fig. 3. The map is drawn by taking the 
longitude east of Greenwich and the latitude of any point as the Cartesian 
co-ordinates of the corresponding point of the map. Fig. 3 is drawn in the same way. 
The defects of the arrangement in fig. 3, considered as representing the shape of 
the continental block, are sufficiently obvious, the chief being the absence of any 
indication of an Arctic ocean, and the almost complete submersion of South America. 
On the other hand, the fact that even tolerable agreement in so many respects is 
obtained from a spherical harmonic analysis of the extremely simple distribution 
detailed in the table of § 56 may be regarded as a confirmation of the theory which 
led us to assume that harmonics of the first, second, and third degrees should be 
predominant. 

Geological Implications of the Theory. 

62. The results appear to admit of a geological interpretation. We have adduced 
dynamical reasons for the hypothesis that the lithosphere consolidated in a shape 
which may be described as an ellipsoid with three unequal axes, with its centre 
of gravity displaced from its centre of figure, and with a wrinkle upon its surface 
expressed by spherical surface harmonics of the third degree ; and we have found that 
the figure of the lithosphere now, as determined by difference of level above or below 

* The coefficients r, c, a, b, c are especially sensitive to changes in the assumed distribution in the Arctic 
and Antarctic regions where the actual distribution is least known. 
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the geoid, is expressible, at least roughly and approximately, by means of harmonics 
of the first, second, and third degrees. Now, if the shape of the lithosphere is at all 
close to that in which it may be presumed to have consolidated, the inference would 
seem to be that, in respect of general features, as distinguished from local irregularities, 
the positions of the continental blocks and oceanic regions have not changed much 
since the date of consolidation. This view has in recent times met with considerable 
support among geologists. 

The theory also enables us to make some attempt to indicate the general nature 
of those changes which could be expected to take place. In estimating the value of 
such an attempt some allowance must be made for the fact that the theory of an 
elastic body in a state of initial stress is very far from complete. We try to follow 
out certain clues drawn from the scanty knowledge we possess of bodies in states of 
initial stress. Amonff these the behaviour of cast iron under tensile tests is perhaps 
important. It is well known that cast iron which has not previously been tested 
exhibits a stress-strain curve which is essentially different from that of mild steel, but 
that, after several tests, its behaviour approaches to that of steel. It has been 
conjectured that the tests have the effect of gradually removing a state of initial 
stress, and thus reducing the substance to a " state of ease." That state of a rotating 
gravitating planet which would correspond to a state of ease in solid bodies at its 
surface would seem to be a state in which the material would be arranged in 
concentric spheroidal layers of equal density, and the external surface would be an 
oblate spheroid, the ellipticity being determined by the speed of rotation and the 
distribution of density ; the state of stress in the planet, when in this state of ease, 
would be one of hydrostatic pressure, and the surface would be an equipotential 
surface under gravity modified by the rotation. The partial reduction of the body to 
the state of ease would be effected by gradual stages, probably of the nature of local 
fractures. Now the wrinkling of the surface, expressed by harmonics of the third 
degree, arose as a consequence of the displacement of the centre of gravity and of the 
ellipsoidal configuration. It would at first be small in comparison with the deviations 
from spherical symmetry which are expressed by harmonics of the first and second 
degrees. We should therefore expect that the tendency of secular change in the 
shape of the lithosphere would be to diminish the coefficients of the harmonics of the 
first and second degrees. An exception must be made in the case of the coefficient e 
of the zonal harmonic of the second degree ; for this coefficient represents a difference 
of ellipticity of the meridians of lithosphere and geoid, and these ellipticities depend 
upon the speed of rotation. When this coefficient is left out of account, the harmonic 
inequality of the second degree represents ellipticity of the equator # and obliquity 
of the principal planes ; the harmonic inequality of the first degree represents 
displacement of the centre of gravity from the centre of figure. If the coefficients of 

* G. H. Darwin concluded from his theory of the tidal deformation of a viscous spheroid that an initial 
ellipticity of the equator would tend to be obliterated. 'Phil. Trans. Eoy. Soc.,' vol. 170 (1879), p. 30. 
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the harmonics of the first degree have ratios anywhere near to those given in the 
table of § 60, the great circle along which the harmonic inequality of the first degree 
vanishes has a pole somewhere in south-eastern Europe and the opposite pole in the 
Pacific Ocean. The inequality is positive in Europe, most of Asia, Africa, North 
America, the northern and central parts of the Atlantic Ocean, and the Arctic 
regions. The effect of a gradual diminution of the coefficients of the harmonics of 
the first degree would be a gradual emptying of the Pacific Ocean, accompanied 
by a rise of sea-level around the shores of the Atlantic Ocean (except towards 
the southern parts of Africa and South America) and around the northern and 
western parts of the Indian Ocean. It has been held that such an effect has 
taken place and constitutes the reason for the difference between a " Pacific coast" 
and an "Atlantic coast," # The ratios of the coefficients of the various harmonics 
of the second degree for the two distributions considered in §§ 56-60 are widely 
divergent, but they agree in leading to negative values for the harmonic inequality 
of the second degree in the regions contained within oval curves which lie within 
the basin of the Pacific, and also in the antipodal regions. In a large part of the 
Pacific region the harmonics of the first and second degrees reinforce each other ; 
in the antipodal region they are antagonistic. Diminution of the coefficients of 
the harmonics of the second degree would be manifested by a fall of sea-level in 
the Pacific, and also in a region antipodal to some part of the Pacific. It may not, 
perhaps, be altogether fanciful to see in the gradual reduction of area of the " Central 
Mediterranean Sea" of Mesozoic and Tertiary times f the effect of a continual 
reduction of those coefficients of harmonic inequalities of the second degree which 
represent ellipticity of the equator and obliquity of the principal planes. Whether 
these conjectures as to the particular regions which may have been affected are 
acceptable or not, it can safely be said that the effects of changes in the harmonic 
inequality of the first degree, and in those of the second degree which we are now 
considering, would be progressive in the same sense at the same place. They would 
be manifested in a tendency of the sea to fall in certain regions and to rise in certain 
complementary regions and gradually to flood wide areas. The gradual character 
of the positive movements of the strand-line, by which wide areas have been sub- 
merged, has been emphasized by Suess. J 

The surface of the lithosphere is nearly an oblate spheroid which does not coincide 
precisely with an equipotential under gravity modified by the rotation; it is less 
oblate than the geoid. The surface of a shallow ocean covering an oblate spheroidal 
planet, whose outer surface is not exactly an equipotential surface, is an oblate 
spheroid, and its ellipticity is a certain multiple of the ellipticity of the surface of the 
planet. The ratio of the two ellipticities depends partly on the rigidity of the planet, 

* E. Suess, < The Face of the Earth' (Translation), vol. 2, Oxford, 1906, p. 553. 
t Mel., pp. 258, 299. 
% Ibid, p. 543. 
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partly on the ratio of the density of the ocean to the mean density of the planet, and 
partly on the angular velocity. Owing to tidal friction, the angular velocity of the 
Earth's rotation is being gradually diminished. The effect of this is that both the 
ellipticity of the lithosphere and that of the geoid are being diminished, and the 
difference of these ellipticities is also being diminished. If, therefore, the shape of 
the lithosphere were continually adjusted to the instantaneous angular velocity, the 
value of the coefficient e of § 57 would diminish continually, and the adjustment 
would involve a continually increasing deformation. Eventually the deformation 
would be so great that the strength of the material would be too small to withstand 
it, and local fractures would take place. # There is, therefore, a constant tendency 
for the sea-level to rise in the polar regions and to fall in the equatorial regions, the 
separation between the regions of rising and falling sea-level being marked by the 
zero-lines of the zonal harmonic of the second degree, that is, by the parallels of 
latitude, about 35° N. and 35° S. This rise and fall would be checked at intervals by 
subsidences, accompanied by series of earthquakes, in equatorial regions. 

The effects produced by diminution of the displacement of the centre of gravity, 
and by changes in the ellipticity of the equator and in the obliquity of the principal 
planes, appear to be of a different character from the effect of diminishing angular 
velocity. The former would seem to be spasmodic and occasional, but always in the 
same sense at the same place ; the latter would appear to consist of continuous 
movements in the same sense, extending over long periods, which are followed by 
comparatively short periods of spasmodic change in the opposite sense. 

These remarks are frankly speculative, and I am well aware that many causes 
which have contributed to geological changes have been left out of sight. They are 
put forward as tentative suggestions which, it is hoped, 'may prove to be of, some 
assistance in the solution of some of the still unsolved problems of geology. 



"y best thanks are due to Professor W. J. Sollas and Dr. H. N. Dickson for 
much kind help in regard to geological and geographical questions. 



* According to a "Note" in * Nature/ vol. 39 (1889), p. 613, this effect of diminishing speed of rotation 
was noted by M. A. Blytt. I have not seen the paper referred to in the " Note," 
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